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Abstract 

In this course, I talk about the source of mathematical constructivism and its role 
in the future development of theoretical physics. I describe what physical construc- 
tivism is and why it is necessary for the penetration of exact methods of theoretical 
physics to the area of complex systems, which formally belong to the others natural 
disciplines. I describe the concrete heuristic for the creating of models of construc- 
tive quantum theory the method of collective behavior. I represent the constructive 
viewpoint on the quantum computer, which treats it as the model object of many 
particle quantum physics, and the practical recommendation concerning its build- 
ing. Due to the known inertia of the educational system constructive methods in 
mathematics remains mostly unknown to the wide physical community, and I hope 
that this course will stimulate the interest to the new possibilities, which these 
methods open. These possibilities represent interest to natural scientists and for 
programmers as well. 
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1 Lection. Hilbert program and its fate 



The strongest source of optimism in the job in science is the perception of its unity. 
This unity has the certain embodiment in the form of the priority of the exact methods 
over the unclear reasoning that we can formahze as the leadership of the mathematical 
methods over all area of natural knowledge. The gradual but continuous progress of 
mathematics to the end of nineteenth century led David Hilbert to his famous program, 
which core was the thesis about the necessity of the axiomatization of natural knowledge, 
and the representation of it in the form of some deductive theory. He expressed this 
aim approximately as follows: At first we put into order mathematics (e.g. prove its 
consistency, and, correspondingly, its absolute reliability). Then we take up physics and 
represent it in the form of mathematical theory. Then we do the same with chemistry, 
and after that with biology. His program contained the various concrete mathematical 
problems, which were solved in the different times, but the most valuable in Hilbert 
program was just its core, which expressed the factual plan of conquest the natural sciences 
by mathematics. 

Hilbert program cracked already at the first stage: it failed to manage with mathe- 
matics itself! In the thirties of nineteen century, Kurt Gedel proved that in any factually 
consistent theory the sentence about its consistency is improvable in it. It means that 
we cannot in principal guarantee that the simple development of mathematics in one day 
will not lead us to the destroying of its main parts in the form of suddenly appeared 
contradiction. The realizing of the real sense of Gedel result took some time but it turned 
many to the shock. It seemed that mathematics has lost the status of fundament and 
nothing could return the perception of reliability, which comes with an exact proof of 
theorem, where each step is clear and causes no doubts. 

For us the most important conclusion from the crack of Hilbert plan is the deep 
understanding of the unity of mathematics with experiments. An experiment represents 
the integral part of mathematics itself, and just it serves as the source of certainty in the 
correctness of results even in the case when these results seem abstract! We can put it 
more exact: physics and mathematics represent the same discipline and their division is 
no more than the question of convenience in the organization of science, at least in the pre 
computer era. This thesis is important for the following and I express it more concrete. 
There is no what somebody call the physical sense independent of mathematics. Just 
the mathematical apparatus determines what has and what has no the physical sense. I 
propose the listeners to think about it themselves, and hope that the following material 
will give the additional arguments for this. 

We thus see that the realization of the plan of axiomatization of natural sciences 
even did not reach the border of the real physics. How can we estimate it? Is any plan 
of joining of the natural sciences the fiction, which does not deserve the attention and 
efforts at all? The main criterion in science differing it from the other forms of mental 
activity, was always the possibility to predict the happening of events. I propose you to 
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compare the level of reliability of predictions in the area with the developed mathematical 
methods, for example, in space mechanics in the scale of Solar system, and the prediction 
of behavior of a living thing, to which mathematical methods are applicable in much lesser 
degree. The difference is so huge that the connection between the degree of applicability of 
mathematics and the reliability of predictions becomes evident. Factually, the all theory 
in the natural sciences is built around the most efficient applications of mathematical 
methods. 

We have no right to refuse from Hilbert program because it would equal to the ca- 
pitulation of the science at all. It is also doubtless that the aim is just in the complex 
system dynamics, e.g. the elucidation of processes, which formally belong to chemistry 
and biology. The question then concludes in the choice of appropriate means for the 
realization of this aim. It makes us to pay the special attention to the mathematical 
apparatus. Can the axiomatic method cope with the functions required from this appara- 
tus? The actuality of this question now follows from that in the modern science deals with 
complex systems where the limitations caused by Gcdcl type results play the key role. 
Factually, the axiomatic methods discredited itself already in the crisis of mathematics we 
mentioned. We then cannot even hope that it can rehabilitate in the applications to the 
real physics; it can serve as the auxiliary tool only, connected with the systematization 
of existing knowledge, but cannot be the working tool of the theory. The perception of 
unreliability of the axiomatic approach leads to the decreasing of enthusiasm in the funda- 
mental substantiation of the observed phenomena among the wide mass of investigators, 
and to the appearance of surrogate specialized theories of this or that phenomenon. In 
any case, the axiomatic method is connected with the abrupt decreasing of the predictive 
force of the theory in its application to the more and more complex systems, so that when 
we approach the borders of biology this force becomes practically zero. 

We then see that the possibilities of the traditional mathematics in the applications 
to natural sciences are close to the exhausting. It brings the idea of the modification 
of mathematical apparatus, for the work with the complex systems. One must clearly 
understand that it does not mean the replacement of mathematics by the kind of philoso- 
phy, or something of that sort, but the modification of mathematics itself, which logically 
follows from its crisis after Gedel theorems. 

In the order to understand the new adventures, which this modification brings, we 
must get acquaintance with the main theses of mathematical constructivism, which we 
intend to use instead of the standard mathematics. 

1.1 Constructive mathematical logic 

Fortunately, the understanding of the dramatic situation in mathematics after the crisis 
led mathematicians to very fruitful turn, which is called constructivism. The principal 
exit from the crisis found in mathematics will soon spread its influence to physics, and 
1 will speak about it. This solution is the constructive mathematics, and we take up it 
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now. Constructivism limits the sphere of apphcabihty of mathematical abstractions by 
the exact computational procedures. It treat as right only such logical formulas of the 
form A OT B for which we can point what is right exactly: A ot B. The sentence A 
or not A is not the universal truth in constructivism. It is the sense of the constructive 
mathematical logic, or intuitionism. Founders of intuitionism were Kronecker and Brower. 
The first of them disputed about intuitionism with Hilbert who could not accept the severe 
limitations that intuitionism imposes to the logical constructions: there are no proofs of 
the pure existence in it. 

Hoverer, already in that times the magic connection between the constructive mathe- 
matical logic and quantum theory became evident. John von Neumann and Garry Birk- 
goff put forward the thesis that the intuitionism is more appropriate for the description 
of quantum theory than the classical mathematical logic. Their arguments rest on that 
the constructive understanding of logical sentences corresponds to the extraction of the 
information about quantum states via its measurement. For example, the formula {A or 
B) in its constructive treatment means that we have performed the measurement and 
know in what state exactly our system is: in A or in B. 

There is very important for physics concept of pluralism in constructive mathemat- 
ical logic that is the plurality of the logical sentences understanding. Models of logical 
intuitionistic theories, in contrast to their classical analogs, contain no absolute true es- 
timations. These estimations may be not only the truth or the false, but also uncertain, 
e.g. not known. In the last case, this estimation can change in time. In the models of 
intuitionistic theories, there is the principal parameter: the time. It puts the construc- 
tive logic in the particular relations with physical experiments, in comparison with the 
classical logic 

The pioneer paper of Neumann and Birkgoff had no big continuation. Theoretical 
physics uses not the logic directly, but the computational instrument, corresponding to 
this logic; in classical case, it is the solution of differential equations, integration, differ- 
entiation, group representations, etc. Any logical interpretation without the obtaining of 
new results has the same secondary value as any interpretation of physical theories at all: 
it only helps to understand already obtained results but gives no new ones. 

1.2 Algorithm theory 

The appearance of constructivism in mathematical logic was immediately stimulated by 
the crack of Hilbert program and the axiomatic method. Nevertheless, in that time the 
main notion of constructivism algorithm yet has not been defined. Turing, Post and 
Markov-young did it independently. Many bodies connect the necessity of the formal 

personally have learnt about the work of Neumann and Birkgoff, written in thirtieths years of 
twentieth on the seminar of A.G.Dragalin. Just Albert Grigorievich explained me the meaning of the 
constructive ideas of A. A. Markov- young whose lectures I listened. However, for the understanding of the 
real value of constructivism I still had to overcome the long way in physics. 
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definition of algoritlims witli tlie need to prove tlie absence of an algorithm with some 
certain property. There is no doubt in the value of such distinguishing with the old 
not algorithmic world. However, the exact definition of algorithms has also the posi- 
tive dimension: it concerns the possibility to compute something practically. Of course, 
many scientists quickly understood the interest for physics (in our country there were Ti- 
honov and Samarsky, the row of their pupils) and the rapid development of what is now 
known as computational methods. These methods master the rich material created in the 
mathematical analysis by the finite differences and methods of linear algebra, and their 
numerous modifications connected with the equations of mathematical physics. In par- 
allel, many specialists developed the notion of algorithms itself (Kolmogorov, Uspensky), 
and elucidate its universality, e.g., the independency from the concrete realization (thesis 
of Turing, Church, Post and Markov- young, which is often called as Church thesis). The 
theory of algorithm complexity began to develop as well (Cook and others). 

However, at this stage the physical applications fit in this narrowly understood sense 
of computational methods. Tihonov and Samarsky ritually fixed this as the known triad: 
model, algorithm, numerical solution (of the differential equation). Of course, this triad 
completely belongs to the classical mathematics, because its leading element the model is 
the classical description of the system by means of differential equations. Algorithms here 
play the role of the technical service. The appearance of computational methods made 
possible to advantage to the borders of applicability of classical mathematics. This area 
includes all classical mechanics and electrodynamics, the thermodynamics, the theory of 
homogenous environments (aero and hydro dynamics, , combined mass processes like 
nuclear reactions in the environment, etc), the theory of gravity and quantum mechanics, 
which in that times deals with quantum effects on the level of one or two particles. 
To nowadays the bulk of principal problems of this sort is already solved but those, 
which nature is purely quantum. We deal just with them. We agree that this area, 
which in nineteenth century belonged to physics (e.g. contained principal questions), 
was completely mastered by the computational methods. The further development of 
algorithm theory seemed predictable but at the end of twentieth century, the surprising 
intrigue appears concerning the quantum computer, which we discuss further. 

The main conclusion, which we must do from the development of algorithm theory and 
its applications, is that all reasonable procedures of theoretical physics are algorithmicJl 

We can treat as simple the problems we spoke about earlier. We can define these 
problems as ones, for which quantum effects, needed for predictions, touch one or two 

^My personal practice shows that the opposite arguments that sometimes appear in the scientific 
literature comes from the poor knowledge of their authors in algorithms. In the detailed analysis of 
concrete situations, which someone treats as non-algorithmic these arguments disappear as a mirage in 
desert. It does not concern the attempts to include to the area of theoretical physics such an object, as a 
human brain. In the dispute about this matter between A. A. Markov- young and A.N. Kolmogorov I take 
the side of the last one. In my opinion, we cannot include this object in the area of physics, at least while 
the constructivism yet does not reach the borders of biology. 



6 



particles. In the other words, these problems are physically simple. 
1.3 Constructive mathematical analysis 

Constructivism represents the fundamental direction. This radically distinguishes con- 
structivism from the computational methods, which bears the applicable sense. Con- 
structive mathematics continued to develop in parallel and independently from the de- 
velopment of computational methods; this development made constructivism more and 
more important for theoretical physics. The next step was the creation of the constructive 
mathematical analysis, which authors were A.A.Markov young, and his pupil G.S.Tseitinjl 

The sense of the constructive mathematical analysis is the refusal from the purely 
classical (and non-physical) notion of a real number. In the reality, we work not with 
real numbers themselves, but with their rational approximations, determined by some 
algorithm. These algorithms we treat as the constructive real numbers. It immediately 
gives us the non-solvability of the identification problem: it is impossible to determine 
algorithmically the equality of two given constructive numbers, though this is not an 
obstacle in the operations with such numbers. All real numbers, which can be determined 
by some reasonable way are constructive, for example, independently of their algebraic 
properties, for example, vr*^, etc. The attempt to find an example of non-constructive 
number immediately requires non-algorithmic processes that have no physical sense. 

We define a constructive function of the constructive real variable as the algorithm, 
which gives us the approximation /„ of the function f{x), given an approximation x„ of 
the variable x. This approach to the definition of functions completely corresponds to 
physics because in the reality we never know the exact values of x or /(x), but only their 
approximations. Of course, if we fix some small threshold for the error e > 0, which we 
agree to treat as admissible, we can use the usual definitions. The principal difference 
appears only if we launch the error to zero. If it is possible to use the actual infinity in 
the classical mathematics, for example in the notion of the limit in the sense of Cauchy, 
in constructivism it is forbidden! Here we have no access to the exact values of x and 
f{x), but their approximations only. It leads to the surprising corollaries. Namely, in the 
constructive mathematical analysis there is the theorem of A.A.Markov the young and 
G.S.Tseitin, which says that 

Any constructive function of a constructive real variable is continuous in 
each point of its domain 

It means the impossibility to use in constructivism non- continuous functions, like 
Heaviside function. Constructivism requires that the complete algorithm describe the 
real situation with all the details with involving all physical laws, which act in the con- 
sidered area of configuration space. For example, considering an electron in the Coulomb 
potential on very small distances we must apply not electromagnetism only but also nu- 

•^The valuable deposit has been done by A.A.Muchnik, Y.Voscowakiks, B.S.Kushner and others. 
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clear interactions. The practical application of constructivism presumes the realization 
of algorithms, which cannot work infinitely. All processes of sequential approximation 
must then finish. For example, we must cut potentials of Coulomb type in some manner 
that we presume implicitly to avoid the divergence of rows in quantum electrodynamics 
that is not a legal way in classical mathematics. In constructivism, this constriction is 
unavoidable, and the reasoning in QED obtains the legal status. 

Moreover, the constructive mathematical analysis completely preserves all analytical 
technique used in theoretical physics that makes legal all its results, e.g., all results of the 
modern physics remains unchanged in constructivism. Constructivism makes legal also a 
reasoning, which needs the complex substantiation in classical mathematics, and some- 
times we cannot even reach such substantiation. These substantiations are not necessary 
in constructivism. Constructive physics does not require the exact values of magnitudes. 
Its aim is the different; we will speak about it further. 

1.3.1 Constructive algebra 

A.A.Markov young started the constructive algebra by proving the theorem about the 
existence of a finite generated semi group with the algorithmically unsolvable problem 
of equality of words. It was the first train of the notion of the normal algorithm he 
introduced (the definition of algorithm equivalent to Turing machines). The connection 
between algebra and algorithms became explicit from this time, through the operations on 
words in the basic algebraic system semi group, and it further developed on groups, rings, 
algebras and the other algebraic systems. □ By a constructive algebraic system, we mean 
a system, which elements are encoded by the natural numbers, all operations are effective 
and the problem of equality of elements is solvable. This definition can be made more 
weak by various ways, for example we can consider the finitely generated systems for which 
the problem of equality not necessary solvable, etc. For the physical applications, it is 
important that matrix algebras, which constructive version uses constructive real numbers 
does not have the solvable equality problem, as constructive real numbers themselves, 
but support the analogy with its classical versions as well as in mathematical analysis. It 
means that all computations with matrix algebras remain valid in the constructive version 
of algebra. 

1.4 Programming and its role in constructivism 

Programming appeared as the art of practical realization of algorithms. We here interest 
in natural problems, and I use the term modeling though the sense we acquire to it is more 
general than usual. By the modeling, we understand the description of the considered 
system evolution in time and the family of such evolution by a computer, such that these 
descriptions we formulate in any terms permitting the computer processing. For example, 

''The further deposit to this direction was done by A.I.Maltsev, P.S.Novikov, S.I.Adian and others. 
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it can be the solution of the system of differential equations; it is the usual understanding 
of a modeling in the narrow sense. It may be also the direct modeling of the process where 
we directly represent its internal microscopic mechanisms as the operations or subroutines 
up to some known physical level (for example, to the level of a single molecule). 

The physical constructivism opens for the programming the new and exciting per- 
spectives, because it obtains the much more high status than in the standard classical 
modeling. The point is that the model of quantum process requires the presence in it 
the specific subject called the observer. Quantum mechanics says nothing about how to 
distinguish this object from the others. We do not know, who or what can and who can- 
not be the observer. We presume only that the observer is the principally more complex 
system than the considered system. Moreover, it is impossible to include the observer into 
consideration as the object, which has some quantum state |^o6s), since it must obey not 
quantum but classical laws. It means that in the entangling with it the state of considered 
object the collapse of the common wave function happens that is the act of observation. 
It means that there is not purely quantum mechanics, because for it something obeying 
classical laws must exist. 

This strange feature of quantum theory always involves the deserved admiration at this 
science. However, the preserving of this feature would make our main aim inaccessible. We 
must have such a model, which contains all the possibilities of scalability to the complex 
systems. It means that we have to include to our model also the mysterious phenomenon 
known as the collapse of the wave function. Our model must work in principle without 
any observer, e.g., the observations (that are measurements) must happen in it naturally, 
in its ordinary evolution. This requirement is very serious and we must explain it in more 
details. 

The uncertainty of the status of observer in quantum theory has the deep source. This 
source is the existence of principally new phenomena, which appears only in complex sys- 
tems but not in simple. The main of such phenomena is the quantum entanglement. The 
other, not known yet phenomena are also possible, which concern the growing complexity 
of the considered system. Just the existence of such phenomena is the reason of division 
of the natural knowledge to the different sciences. We must take into account it in the 
realization of our project. This is why we cannot ignore the traditional structure of quan- 
tum theory. However, it turns that in the mathematical constructivism there is already 
the certain place for the observer! This is the place of the oracle in computations. 

Let we be given an integer function F{x) of the unknown nature (for example, it can 
be non-computable, e.g., have no algorithm computing it). What means the computation 
with the oracle F? We imagine that the algorithm, accordingly which we perform com- 
putations, can in some instants query the oracle with the question: what is the value of 
the function F on the number x, which contains in the certain register of the memory? 
After query, the algorithm temporarily stops the computation and waits the reply from 
the oracle. After the receiving of the reply, it puts it to the other certain register and 
continues the computation. Then it can repeat the query on the other word x etc. It is 
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the computation with the oracle. Thinking a Httle, we guess that the user of the program 
plays the role of the oracle. We then conclude that the user of the simulating program 
plays the role of the observer of the quantum system. The roles in the modeling are thus 
completely determined and in principal, we can begin to build programs. 

We must make the important notion about the role of the simulating programs in 
constructivism. The known fact is that it is impossible to predict the work of algorithm 
even on one-step ahead. There is the single one way to learn what an algorithm does: 
to launch it and observe its work step by step. It radically distinguishes an algorithm 
from its particular formula, where it is typically possible to predict the behavior 

of the function it determines even to the infinity. This difference radically changes the 
criterion of success of a theory and the form of a theory in constructivism. We cannot have 
as a criterion the exact computation of some physical magnitude found in experiments 
that is the tradition of the usual physics. The criterion of validity of a physical theory in 
constructivism is the rightness of the dynamical scenario built by the modeling algorithm. 

A physical theory in constructivism is thus a heuristic for the building simulating algo- 
rithm, and the main object is the dynamical scenario. A theory is valid if this algorithm 
gives the right video film of the considered process. This is the success criterion in con- 
structive physics. The exactness of the estimations of separate magnitudes is important 
only it helps to build the right algorithm, but not more. The exactness of estimations 
plays then the secondary role. 

The new role of programming is evident from this. In constructivism, it is the work- 
ing tool of the theory. Its role is the same as the role of analytical computations in 
the traditional quantum physics. It brings the new requirements to programs, which 
constructivism development requires. These programs must be universal, not narrowly 
specialized. They must allow the simple introduction of the new particles and interac- 
tions in the considered system and permit to work simultaneously the different groups of 
specialists. Such programming seriously differs from its traditional treatment as the tech- 
nology of the realization of complete mathematical solutions, and it requires the special 
attention. 

1.5 What systems are complex 

We now explain what we will mean by a complex system. It is impossible to predict 
the evolution of such system by systems of differential equations or something of this 
sort. This definition is not satisfactory because it is non-physical. The impossibility 
to describe by systems of differential equations appears as the random character of the 
behavior of the system, even in the case when its usual behavior admits such description 
but there arc singular points in which the evolution acquires the random character. For 
example, the behavior of the media can be in the critical dependence of a few molecules or 
atoms. Here we distinguish two cases. In the first case to get the complete description, we 
only need to apply classical mechanics to these atoms. In the second case, the complete 
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description requires the application of quantum mechanics. It happens if the mechanical 
movements in the macroscopic volumes depend on chemical processes that involve a few 
molecules. We have no experimental evidences that these two types of complex systems 
are independent from each other. Moreover, we can certainly state the opposite. Namely, 
the source of any randomness in the Nature is some quantum phenomenon. 

One could object, on the example of the collisions of billiard balls, which in the case 
of large numbers of collisions give the completely chaotic picture, whereas the collisions 
are classical. The point is that the classical description is effective for a few sequential 
collisions only. If the sequential collisions are more than ten (and the number of the 
degrees of freedom is more than one otherwise there is no place for the randomness), the 
classical description leads to the complete chaos. The point is that the error in targeting 
grows exponentially with the time that requires the corresponding increasing of accuracy, 
and it quickly leads to the necessity to consider the separate molecules on the surfaces of 
balls and their bounds inside of the ball, and the problems goes out of the framework of 
classical physics. 

We conclude that the systems representing the biggest difficulties for the simulation are 
quantum systems. The computational complexity then always has the physical basement. 
It is impossible to overcome quantum nature of complex systems. The theory of complex 
systems is thus the theory of quantum systems with many particles. 

1.6 Why physics needs the constructivism 

We thus come to the main question about the connection of constructivism with physics. 
We formulate the conclusion definitely: to spread theoretical physics to complex systems 
physics must become constructive. This statement of question concerns the great project 
of quantum computer and it has the certain sense. 

1.6.1 How quantum computer appeared 

We can really call the project of quantum computer (QC) the great because it decisively 
violates the usual order in the relation between physics and mathematics. In order to 
understand how this project lead us to the necessity of the constructive modification of 
quantum theory we need its detailed consideration. The source of QC is the same idea of 
Hilbert about the mathematization of the natural knowledge, and about the spreading of 
mathematical methods to complex systems. Since the way of axiomatization is closed, we 
have only one way: the simulation, or the modeling. The question is only what algorithms 
are appropriate for this. The building of the simulating algorithms is very serious thing, 
which we call the heuristic. The mathematical apparatus determines it. The traditional 
apparatus of quantum mechanics is classical mathematics, namely, the theory of Hilbert 
spaces. Here the state vector in the space of the dimensionahty exp{m n) determines the 
state of the system, where n is the number of real particles in our system, m is the number 
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of grains in the configuration space accessible for the simulation, so that m — 1/g, g is the 
spatial resolution grain. The operator of evolution without measurements represents the 
unitary operator Ut, which depends on the time in this space. To simulate the evolution 
in this mathematical terms we need to work somehow with the matrix of exponential 
dimensionality that is impossible yet for three quantum particles (the case of two particles 
is peculiar; in quantum mechanics it is reducible to the case of two independent particles, 
as well as in classical mechanics). 

The first who recognized the fundamental obstacle was R.Feynman. He proposed to 
use for the simulation not the classical but the quantum computer. Quantum computer is 
the binary device on quantum bits, operations on which imitate the evolution of the real 
quantum system. This was the revolutionary idea violating the formed status-quo between 
physics and mathematics. At first, it turns that the notion of fast algorithms depends on 
physics. At second, algorithms themselves obtain the new power over physics, and it will 
never be the former science. The last thesis appeared not immediately, but after about 
ten years of rapid development of the mathematical quantum algorithm theory. 

At the first stage of the development of the quantum computing, it became clear 
that Feynman idea is blameless from the viewpoint of the standard formalism, which he 
possessed. Really, Wiesner and Zalka showed that it is possible to simulate the evolution of 
the abstract state vector in the space of exponential dimension by the quantum computer, 
factually operating with the quantity of qubits proportional to the volume of the system 
at hand. The time of such simulation we can make as close as desired to the real time 
of the modeled system that requires the simple change of the constant. Moreover, it 
turned that the quantum computer is able to solve also purely mathematical problems, 
for example to find the expansion to prime integers in the polynomial time (Shor) and 
the search problem in the time \/iV where is the time of the classical search (Grover). 
Many other results appeared then formed the ideology of abstract quantum computing. 
For example, it turned that the bulk of classical computational tasks cannot be sped 
up on quantum computer at all. This underlines the peculiar character of quantum 
computations. There are not the usual computations in the classical sense. There are the 
special information procedures, which in some cases make possible to obtain the result 
much faster than any computation on a classical computer. Just this feature violates the 
usual non-formal bounders between physics and mathematics. Quantum computational 
procedures factually use some peculiar resource, named the entanglement, which does not 
exist in classical physics and in classical computations. 

1.7 Quantum computer as the model object of the new physics 

If a quantum computer was built in its scalable version (versions with a few qubits already 
work, but this is the big difference), this course would never appear. The experimental 
work at its creating was the serious and sufficiently long that we can make some conclu- 
sions. Here the success would mean that the entanglement is the unlimited resource, from 
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which we can extract the information factually without the deep analysis of its sense. 
The algorithm of Zalka and Wiesner makes possible to simulate complex chemical and 
biological systems and use all the advantages that it gives, without the deep analysis of 
the nature of the phenomena accompanying their evolutions. 

However, the different has happened. Quantum computer project slipped and buried 
hopes to conquer complex systems by means of the standard quantum theory! Formally 
speaking, quantum theory has not pass the test for effectiveness in the area of many body 
systems because experiments show that the role of decoherence substantially exceeds the 
effect of abstract quantum computing. For complex systems, the contra intuitive and non- 
mathematical part of quantum theory turns more valuable than the exact mathematical 
part, connected with quantum computations. We could predict this. One of the first who 
understood it was K.A.Valiev, who started the work in the quantum computer physics 
in our country. What conclusion can we make from this? Could we doubt in quantum 
theory itself, despite of that just it led us to the understanding of the real problems of 
the Nature? Wc could not do it even if it would not such a triumph of this science in 
the structure of atoms and molecules, optics and all the other happened in physics in 20 
century. To make the right conclusion we should return to the thesis about the leading 
role of the mathematical formalism, which we discussed earlier. The conclusion can be 
only one. The classical mathematics was wrecked. Its pretentions to the conquest of the 
world of complex systems turn vain. The applications of algorithms in the traditional 
methods, giving the effect for classical systems turned out too week in the principally new 
situation when we study quantum physics of many body systems. 

All physical applications of the theory of differential equations have the certain sphere 
of applicability connected with the existence of the elementary quantity (grain) of each 
magnitude. The derivations of differential equations (oscillations, heat transfer, etc.) 
themselves rest on the application of the main law of interaction, which is limited by 
this grain. If / is a quantity with some physical sense, it has the minimal nonzero value 
A/. We consider quantum theory from this viewpoint. The wave function |^) of any 
system must be then grained, e.g. it has the grain 5$. We now consider the wave function 

N-l 

of n particles of the form l^') = ^ Aj|j) where N — 0{exp{n)) is the total number 

of basic states of the considered system. The typical value of |Ajp then equals N~'^. 
The detection of this magnitude requires the time of the order N , and it is impossible 
to estimate it practically already for the relatively simple systems of a few particles, let 
alone the practically interesting case of bio molecules. It means that we cannot in principal 
guarantee, that quantum states we create are really that quantum theory predicts. Let 
us suppose that the nonzero value 5* really exist and is sufficiently small to make no 
influence to quantum theory predictions for one or two particles (the case of two particles 
is reducible to the case of one particle). If we then increase the consider system, the 
typical value of amplitudes decreases exponentially, and we reach the threshold 5^^ very 
soon. The behavior of slightly more complex system will then differ substantially from the 
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solution of Shredinger equation, which the experimenter treats as decoherence. We can 
show (see [1]) that the existence of the amphtude grain 6^, immediately gives Born rule 
for the computing of quantum probability. This is the constructive approach to quantum 
theory. Its specifications give us the mechanisms of modeling of the evolutions, which we 
call constructive heuristics. 

How can we help to build a quantum computer? We must study its physics. The 
physics of quantum computer is the key to its practical building. We do not know to 
what extend we can scale it, e.g., how much is the value 5* in the reality. Quantum 
computing seems now possible up to few qubits, may be up to tens qubits. However, 
where the border lies behind which decoherence makes the scalability impossible? I have 
not answer to this question. Do we need to develop quantum computations? Certainly, 
we need. We must understand how to realize basic principle of quantum computing on its 
real models with decoherence. It helps us to recognize its real borders and its nature; in 
the other words, we then hope to understand where the new physics of complex systems 
begins, which I spoke about earlier. 

However, we now cannot limit our efforts by the movement in one direction only. 
Classical mathematics discredited itself just in the area, where the support of mathematics 
is most important: in the theory of complex systems. There is no exact description of 
decoherence in the standard apparatus, whereas this phenomenon plays the key role in the 
physics of quantum computer. What solution of this problem proposes the constructivism? 
The constructive approach binds us to build algorithmic models of quantum evolutions on 
classical computers. Decoherence we must then treat as the unavoidable constriction of 
complex quantum states that happens due to the deficiency of the classical memory, which 
rapidly increases when the considered system growsjf] The building of such algorithms 
and programs requires the special heuristic, which we discuss further. 

We thus have to rest on the constructive mathematics. This presumes the serious 
change of the standard criterions and views. The main in this unavoidable matter is to 
get into the way of pluralism and take it as the native property of things. This side of 
constructivism is determinant for the success. Further, we see how to make quantum 
theory constructive. 



^Of course, the rest on classical computers does not mean the refusal from the usage of the limited 
quantum processors, for example, for the simulation of entanglement. 
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2 Lecture. History of collective behavior 



We saw that the main criterion in constructivism is the rightness of the video film about 
the considered process. This criterion is vahd only in the case when the algorithm cre- 
ating the video film is the general, e.g., it can process not the single but the wide range 
of systems. The modeling programs then must be universal. The universality in the pro- 
gramming requires the using of the instrument for developers that is the programming 
complex, which makes possible to elaborate the simulation programs for the particular 
cases. Speaking about the simulating algorithm, we will thus mean just this instrument 
for developers. Our algorithm must rest on the universal methods and notions and its 
structure must be completely transparent. The user must have the possibility to change 
the input data and check its work, etc.; we must guarantee the complete verifiability of 
the results. The structure of algorithm must be then universal and rest on the certain 
modification of the standard quantum theory. 

This modification we will describe in the next lecture. It is not the interpretation 
of quantum mechanics. It is the constructive truncation of quantum theory. It means 
that we reduce the formal capabilities of the Copenhagen quantum mechanics that makes 
it coordinated with the constructivism. We call such modification the heuristic. The 
different heuristics are possible depending on what lays in its basement. For example, 
one can formulate the heuristic so that it applies classical mechanics everywhere but 
the processes, which involve electrons, and for them it apphes the wave functions. This 
way could give the initial advance to the chemistry, because this approach is typical in 
quantum chemistry. Nevertheless, this simple approach is not appropriate for our aims 
because it does not satisfy the property of universality. 

The method of collective behavior rests on the representation of one quantum particle 
in the form of the ensemble of point wise particles. This is not the new idea. The most 
known and successful allied approach is the formulation of quantum mechanics in the 
form of Feynman path integrals. 

2.1 Feynman path integrals 

The remarkable principle of quantum mechanics is the principle of correspondence, ac- 
cording to which any physical magnitude corresponds the hermitian operator acting in 
the space of states of the system, where this correspondence ensures the passage from the 
quantum mechanics to classical for the value of action large in comparison with Plank 
constant. There is no quantum physics without classical in the strong sense of the word. 
There is not purely quantum theory; it must contain classical physics because the mea- 
surement procedure is impossible without the measuring device obeying classical laws. 

This impossibility to separate quantum physics from classical has the other sense 
concerning heuristics. The heuristic in the standard quantum theory as well as in the 
classical physics is non-formalized system of notions and agreements, which determines 



15 



how we must apply laws and formulas to obtain the result having the physical sensqj. 
The classical heuristic stands behind all advantages of quantum theory including the 
electrodynamics. This heuristic rests on the notion of a point wise particle. 

The simplest form of systems is those contain only one particle or is reducible to one 
particle in the sense that the approximation of their dynamics by the simple combination 
of one-particle systems is satisfactory. The example is a particle in some potential. A 
real potential is the sum of deposits from quanta of this potential - photons. Speaking 
about a particle in the potential, we mean the approximation resting on the following 
agreements: a) there is the large total number of photons, and b) the emission of one 
photon does not cause the change of the source. We consider two interacting particles 
with coordinates ri and r2. If we ignore quanta of the field again, we could introduce the 
new coordinates R = (ri + r2)/2, r = ri — r2 that reduces our system to two independent 
particles, e.g., to the simple combination of one particle systems. For three particles, this 
trick is impossible and the case of three particles represents the kind of a touchstone for 
the checking of various hypotheses in quantum informatics. 

In this paragraph, we consider systems reducible to one particle. For such systems in 
quantum case the classical heuristic is valid, which allows the reformulation of all conclu- 
sions of quantum physics on the language close to the classical (we call this language quasi 
classical). This heuristic allows the using of such terms as trajectory", "the movement of 
particles along the trajectory", "deposit of the trajectory" etc., despite of that the formal- 
ism of wave functions contains no trajectories and no deposits. This powerful heuristic 
stands behind all advantages of quantum theory achieved to nowadays. Moreover, the 
constructive physics requires this heuristic because we have no alternative. R.Feynman 
did the first step on this way, and formulated quantum mechanics on the language of the 
so-called path integrals. 

The idea of this language is easy. Let us consider the flight of the particle from the 
point 1 to the point 2. It can fulfill this flight along the different trajectories. We accept 
that there is some algorithm determining for a given path 7 the number ^(7), and for 
the different such numbers Ai, A2, . . . ,Ah their result K{Ai, A2, . . . , Ah), which is in turn 
the number. The more the module \K\ is, the more probable our particle will occur in 
the point 2 provided it was initially in the point 1. This idea contains no numbers but 
only evaluations. However, it points to the algorithmic scheme answering roughly to the 
question where will pass the particle from the point 1. We conclude that a quantum 
particle factually represents not one point travelling in the space, but the set of points, 
each of which has its own path. These points are equitable in the sense that each of them 

^The analogy with a juridical system is pertinent where besides laws there are the procedures of their 
applications. It can be sub legislative acts, or precedents, or something else. It is important that this 
addition is unavoidable because laws will not work without it. 

^We are always speaking about the electromagnetic potential. However, all we say about the quantizing 
is right for the other potentials including the nuclear and gravitation. For the last, the question is only 
how to divide it correctly to quanta. 
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equally represents the initial particle. It gives us the right to call them samples of the 
particle. 

A quantum particle looks like a galaxy, which stars are its samples, though this analogy 
requires some efforts in the description of interactions. 

It influence to the mixed consideration of systems by Born - Oppenheimer method. It 
considers atomic nuclei as classical, and electrons as quantum objects. 

In view of the above mentioned I call such a set of samples of the same real quantum 
particle a swarm that underlines its unusual character, where each sample represents the 
whole particle. It brings the question: what forces samples to hold together or what 
happens if some sample flies to far from the main part of the swarm. The answer pre- 
sumes that we treat samples as non-erasable and the existence of some mechanism of the 
returning to the swarm distant samples. We consider this question further. 

It would be logically to treat each sample as non- erasable: that is to ascribe to each 
sample its own history. We return to this thesis further. Now we only treat samples as 
the auxiliary tool for the description of the wave function \E', and will redefine them anew 
in the short time frame 6t, on the basement of the wave function "^{t). This swarm we call 
wave swarm because its samples have only the history limited by the duration 
St. Our scheme then looks as the iteration of three main steps: 

• computing of the wave function from the states of all samples of the wave swarm, 

• new definition of samples in the wave swarm, 

• free flight of samples and change of their amplitudes. 

What is exactly the probability of the passage 1 — > 2 ? This question concerns the 
interaction between samples of the swarm. The right answer is known: the probability 
equals the squared module of the value K, and we will estabhsh why it is so. We are 
not bound by the necessity to consider the history of each sample and treat them as non- 
erasable; we accept that each sample is the carrier of some special number corresponding 
to it, and call this number its amplitude. To specify the details of this scheme we have to 
make more exact the meaning of the terms K{Ai, A2, . . . , A^) and ^(7). The first term 
we define simply as the sum 



h 



K{A,,A, 



The second as 



A(7) = ^exp(l%]) 



where A is some constant and 5" [7] is the action along the path 7, defined as 
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where Lagrangian of the considered system is L = E^in — Epot the difference between 
kinetic and potential energy. 

The function K is thus complex valued and we can express it by the formula: 

K{2A) = -! exp(^5[7])/^7 (1) 

^ JT(2,1) 

where the integration on means the summing on the set of all paths T(2, 1), going 
from the point 1 to the point 2. The wave function in the moment ^2 we can express 
through the wave function in the moment ti by the formula 

^(3:2,^2)=/ K{x2M\Xi,ti)-^{xi,ti)dxi (2) 

J Xl 

for any ti, t2, a^2- Here the point 1 has the coordinates xi, ti, the point 2 - the coordinates 
X2,t2. By the way, the formula [2] reflects the fact that the wave function is the value, 
which squared module is the probability density to find the particle in this point. This 
value equals K provided the particle initially was in the point 1 (in the right side stands 
the action of delta- function on the wave function that immediately gives the desired). 
The formula [1] formalizes the rule of computation of K, which we gave earlier. 

Resuming the previously mentioned, we can define all three items in the evolution of 
the wave swarm. The computation of the wave function in the given point is the summing 
of amplitudes of all samples occurring in some cube with the center in this point. The 
redefinition of samples goes as follows. We divide the value of the wave function in the 
given cube to many peaces, e.g., represent it as na, where the natural n is large, and 
create n new samples. We ascribe to each of these samples the random speed from the 
uniform distribution on the large cube. At last, the free flight of samples goes accordingly 

iAS 

to Galileo law, when all amplitudes are multiplied to e~~h~, where AS is the action of 
the sample along the straight line when the time 6t of its flight is fixed. To simplify the 
computations we agree that a is the same for all points and n = n{x) proportional to 
in the given point x. 

The formalism of path integrals is equivalent to the standard. For this we follow the [3] , 
and find the value of the wave function defined accordingly to [21 in the next time instant. 
We need the explicit representation of the kernel K. We suppose that for small values of 
the period e there is only one straight path of the form 1 — > 2, and the integration in[T] 
turns to one summand only. This trick requires the substantiation from the view point of 
standard mathematics, whereas in the constructivism it is legal, because if we limit the 
grain of spatial resolution by some threshold all trajectories become broken lines, which 
sections are straight lines and we obtain the desired. Analogously, it is possible to assume 
that the action equal the product of e and Lagrangian taken in some intermediate point. 
We then have 

^{x,t + e) = jl expC-^L{^,^))^{y,t)dy (3) 

R 
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Substitutions L = mxl/2 — V{x^ t) and y = x + 5 give 



^{x,t + e)= I ^exp{^)exp{--V{x + 6/2,t))^{x + 6,t)d6 (4) 



R 



We see that the main deposit comes from the 5 of the order ^J ehjm. If we expand \I/ to 
degrees of e, keeping the summands of the first order, we must keep the summand of the 
second order of 5. It gives with this accuracy 

R^ 

and in our approximation with the equation 



e^-^V2.e^5 = f ^ (6) 



R 

we have 



^=f?^^V^ (7) 



Now we apply the known integral 



A 

A m 



R 

and obtain 



^ ^ d"^ _ ^ ze^^ he 

dt h 2im dx"^ 

which immediately gives Shredinger equation. 

The method of path integrals thus represents the version of quantum mechanics. Its 
practical implementation presumes the ingenious computations with integrals and algebra, 
and it gives the new insight into quantum theory in comparison with Shredinger equation. 
For example, it is possible to find the kernel of particle in the different fields (see p]), for 
the free particle, it has the form 

K{x, t, 0, 0) = f ^1 e^™^"/2^* (10) 
\ m J 

which shows the character of the movement of its samples to the points x with the speeds 
V = x/t that corresponds to the classical picture of the flight of pieces after the "explosion" 
of the particle initially located in the center of the coordinate system. 
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Path integrals allow the deduction of the uncertainty relation of the type " coordinate 
- impulse". It follows from the consideration of the passage of the particle trough the 
narrow slit of the weight 2b. Slightly clumsy calculations of the kernel of the passed 
particle (see [3], pages 63-64) show that after the passage through the slit the support 
of the wave function obtains the widening ht/mb, which means the obtaining of the 
additional uncertainty in the impulse as h/b, that gives the uncertainty principle of Born- 
Heisenberg: 

5pSx = 2h. (11) 

In the formalism of path integrals the notion of sample plays the secondary role since 
each sample of the real particle preserves its history in the short time frame 5t onl}|§. 
The main role belongs to the wave function, which factually determines how many and 
what samples we must create in each spatial cube. However, just the free movement of 
samples is the step of evolution creating the unitary dynamics. The secondary character 
of a sample in the wave swarm follows from that the density of such swarm is proportional 
to the module of the wave function, not to its square. In the other words, the density 
here if not the density of probability of finding the particle in this point. The brevity 
of the sample history leads to the absence of the dynamical characteristic of the sample 
themselves. The mass of sample, which in this approach is equal to the mass of the real 
particle occurs only in the evolution of its amplitude through the action. Despite of the 
ephemerality of the wave swarm samples their application gives much. In particular, we 
can by means of path integrals formulate the criterion for what dynamics must we apply 
for this particle: quantum or classical. 

Indeed, let us consider the passage from the point 1 to the point 2 along two paths: 
7c/ the path, which is the solution of the classical dynamics equations and ■jnonci, which 
is some other path. Without loss in generality we can accept that the samples preserve 
their history travelling from 1 to 2, e.g., that these point are sufficiently close to each 
other. We compare two deposits to the wave function: 1) the deposit from amplitudes 
of samples flying along trajectories close to 7c/, and the deposit of samples flying along 
paths close to 'jnonci- We denote them by Kd and Knonci correspondingly. Let Sd be the 
action along the classical path. We suppose that 

|5c/|>/i. (12) 

We can suppose that the change of action on the order of value equals to the action 
itself if the path have the general form. Since the classical path 7c/ is the solution of the 
equation = (the principle of least action) the change of the phase on all paths 
close to 7c/ is small. The change of the phase on paths close to "inonci is large because 
of the inequality [T21 because here the equation ^^j^ = does not take place. It means 
that the deposit on absolute value is much more than the deposit Knonci, due to that 
in the first deposit we sum the values with the approximately equal phases, whereas in 

®There is no explicit notion of sample in the original text [3| at all. 
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the second the phases are different. We thus obtain that all the amplitude concentrates 
along the classical path of the particle, e.g., it behaves as a classical particle. Whereas if 
[T2] is not true, the situation changes, and the deposit Knonci can compete with Kd, that 
is the real particle will not move only along the classical path, and will reveal quantum 
properties, for example, interfere in the passage through two slits, etc. We then can accept 
that the inequality [12] is the criterion of the applicability of classical mechanics. We also 
note that the smallness of the action can follow from the lightweight of the particle, small 
speed, or small period. Even massive particles moving slowly in the short time frames 
reveal quantum properties. In the practical description the period St we take such that 
it gives the pithy general picture (not only the isolated particle dynamics). Hence, for 
example, electrons typically reveals not as point wise objects, but as the wave functions, 
whereas nuclei are classical objects. This is Born -Oppenheimer model. It is convenient 
for the cases as atomic physics, where the subject is many electron states, or molecular 
dynamics, where the stretching and rolling of molecules are in focus. In opposite, this 
model is not always appropriate for chemical reactions where the quantum character of 
the movements of nuclei plays the key role. In the same reaction often is necessary to have 
the classical and the quantum consideration of reagents. Born -Oppenheimer model is 
convenient mainly due to its simplicity resting on the large (in about 2000 times difference 
between the masses of electrons and the proton. 

2.2 Monte Carlo method 

For the search of the ground state (the eigen state with the minimal energy) the most 
convenient is the method of Monte Carlo, which we can call the method of static diffusion. 
This method represents a quantum particle as the swarm of its samples so that its density 
equals to the module of the wave function: p = |\E'o|- It shows the narrowness of this 
method because it cannot serve as the probability model of quantum dynamics where 
the density must equal the square of the wave function module. The deep source of this 
difference is the static character of Monte Carlo method, which is not appropriate for the 
modeling of dynamics. The dynamics requires the description of excited states, not only 
the ground state. However, Monte Carlo method is simple in usage and we can treat it 
as the good starting point for the creating of the real dynamical model. 

We consider the following process of the swarm evolution. Each sample with some 
small probability p jumps to Ax along one of coordinate axes to the positive or negative 
direction with the equal probability p. With the probability 1 — 6p a sample stands at 
place. Let we be given a scalar function V on the configuration space. We also accept the 
following rule of creation and annihilation of samples. Let the following process happen 
for each sample with the probability proportional to V. If < 0, then this sample 
generates the new sample located in the same point, if > 0, then this sample eliminates 
some other sample located on the distance lesser than Ax. We call this process the static 
diffusion. Samples have no speeds here and only their coordinates participate. 
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Now we make in Shredinger equation 



9* 

ih— = -—A* + 
dt 2m 

the formal replacement of the variables t — —it. The equation then acquires the form of 
diffusion equation: 

h— = A* - 

dt 2m 

which just describe the stationary diffusion process. If we put the energy levels into 
order of the growing energy values then the evolution of the state vector expanded to 
eigenvectors 0„ will be: 

n 

The expansion of the diffusion equation then obtains the form 

n 

We see that for t — > oo the diffusion process converges to the ground state, because the 
rapidly decreasing exponentials suppress all other states. It is known that the ground 
state contains no phase differences, e.g., it coincides with its module. Hence, to find it, 
we have to determine the initial swarm density distribution and launch the stationary 
diffusion process. It stabilizes on the distribution proportional to the ground state. Of 
course, we should tae care of the conservation of the total number of samples in some 
reasonable frameworks. We can easily guarantee it by the addition or elimination samples 
uniformly accordingly to the existing density. 

DMC method can be easily generalized to the case of n particles. A sample here will 
be a cortege of n samples, each from the swarm of the corresponding real particle, and 
the configuration space will be i?^". 

We note again that DMC method is aimed to the search of the ground state only, 
e.g., for the stationary modeling. To find the dynamical picture we need the method of 
dynamic diffusion, which we describe further. 



2.3 Bohm approach 

The method of Bohm uses the notion of pseudo potential. We identify the square of the 
wave function module with the density of the swarm and its phase with the classical 
action. We then have ^ = p^/^e"^/'*, and 1/m grad (t){f) we can regard as the density 
of the swarm stream. Shredinger equation then becomes equivalent to such system of 
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equations]^ 

^ + div {p/m grad (p) = 0, 
f+ ^{grad<P)^ + V + V^ =0, 

where the quantum pseudo potential Vi = ^{Ap/p + {grad p)"^ / p^) depends on the 
density of samples with the singularity in the initial point. These equations coincide with 
the equations of classical particles dynamics if we accept that Vi has the form of some 
potential. 

2.4 The drawbacks of analytic approach to collective behavior 

We sum up the tricks using the idea of collective behavior: the exact reformulation of 
quantum mechanics on the language of Feynman path integrals, Bohm approach and 
diffusion Monte Carlo method. They are different in aims but have the common idea of 
the representation of a quantum particle as the swarm of point wise particles. They all 
rest on the classical mechanics that determines the limit of their applicability. 

We begin to analyze them with the most valuable Feynman path integrals. In the proof 
of equivalence of this method with Shredinger equation, we use the explicit representation 
of trajectories in the form of straight-line segments on the small distances. It is natural 
for the constructive version of collective behavior, which determines the dynamics of the 
separate samples by the sequence of such steps: 

• free flight of all samples in the period At, 

• interaction of close samples resulting in the change of their speeds by some local 
rule. 

This is the steps of constructive method of collective behavior, which we call the method 
of dynamic diffusion; we describe it in the next lecture. However, from the view point 
of classical mathematical analysis the supposition about the straight form of trajectories 
on the small periods cannot be proven and we must take it as an axiom. Moreover, the 
analysis of quantum trajectories by path integrals (see [3]) shows that the main deposit to 
the amplitude belongs to nOn-regular trajectories, that is the paths, which have the more 
differences in speed the less spatial grain we use for their representation. This property of 
the fast change of speed on quantum trajectories we can illustrate using the uncertainty 
relation coordinate impulse: the less is the grain of spatial resolution, the bigger is the 
dispersion of speeds inside of each grain. It cause difficulty in the formal substantiation 
of path integral method in the classical mathematical analysis!^ 



^These equations derived Madelung. Bohm lied them in the basis of his interpretation of quantum 
theory; we can treat it, with Feynman path integrals the pre image of the collective behavior method, 
which we study further. 

^"Difficulties in the substantiation of formal methods, which give the maximal effect in theoretical 
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Diffusion Monte Carlo method gives good approximations for wave functions of ground 
states. However, it is not applicable for the description of dynamics and is not thus 
completely constructive. Samples in DMC method have no speeds, only positions and it 
determines limitations of this method. 

Bohm approach contains the serious drawback from the constructive viewpoint. The 
mechanism of interaction between samples, which could give this potential, is unclear. 
Moreover, the singularity in the initial point means that in the areas with the small 
density some huge force acts on samples with the unpredictable direction that makes 
impossible to apply this approach to the computer simulation of quantum dynamics. ["1 

The deep source of this difficulty lies in the statistical nature of the wave function. Its 
experimental finding always requires the large number of the repeated experiments and 
the linearity can reveal exclusively as the limit result of them if their total number goes 
to infinity. The refusal from the simple representation of linearity is not thus the serious 
drawback. Moreover, the effective modeling requires the refusal of something more. We 
need the simple mechanism of the interaction between samples in the dynamic diffusion 
swarm. It turns that for the obtaining of such mechanism we must refuse from the 
relatively beautiful systems of differential equations, as was written earlier. It is possible 
to write the system of differential equations only for the fixed spatial resolution. This is 
the price, which we must pay for the method of building of the effective algorithms, e.g., 
this is the price of the constructive modification of quantum theory! 

At last, the important deficiency of classical versions of the collective behavior that we 
have studied here is the lack of the individual history of samples. Factually, it is true also 
for the DMC method because we have to renormalize the wave function periodically, which 
means the addition or annihilation of samples. We further return to the discussion of the 
role of individual history of samples. Now we only note that the absence of the individual 
history deprives classical versions of the collective behavior the additional predicting force 
in the comparison with Shredinger equation itself. Using of the classical mathematical 
analysis as the main tool in theory makes impossible to us to go out of the framework 
of standard problems contained in traditional courses of quantum mechanics, let alone to 

physics are typical. We can recollect generalized functions, asymptotic rows, etc. The typical reaction of 
physicists to these difficulties is the following: let mathematicians take it up; it is their job. This position 
is pragmatic but narrow sighted. The effectiveness of the usage of mathematical apparatus depends on 
the accuracy in the substantiation of used methods, as the secure driving depends on the satisfaction 
of traffic rules. The violation of these rules leads to the crashing results immediately when the usual 
situation: for example, in the area of many particle systems, containing quantum computers. 

^^D.I.Blohintsev treated in [1] the impossibility of the simple interpretation of linear property of 
Shredinger equation solutions in hydrodynamic terms as the main drawback of this way. I think that it 
is not a lack at all. It is impossible to overcome the contradiction between the stochastic nature of the 
wave function and linear properties of Shredinger equation in any ensemble interpretation of the wave 
function. However, it is not required for the building of effective simulating algorithms. Any method of 
matrix algebra has the genetic lack in the usage of computational resources; for the simulation of many 
particles, we unavoidable have to choose between some version of ensemble representation and quantum 
computer. This is why the sacrifice of the algebraic beauty would be reasonable and necessary here. 



24 



study complex systems. 
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3 Lecture. Dynamic diffusion swarm 



Wc now describe completely constructive form of collective behavior, which can serve as 
the heuristic of the modeling algorithms. 

We call it the method of dynamical diffusion swarm. It differs from the previous 
attempts in that we set the aim not the obtaining of the differential equation but the 
mechanism of interaction giving quantum dynamics by the shortest way for the realization 
on a classical computer. 

Why the dynamical diffusion swarm is better than the solution of Shredingcr equation 
by one or the other method? In the solution of a differential equation, we factually use 
Riemann scheme of integration. Computations of the wave function go in all configuration 
space independently of how constructive the interference of amplitudes is. On the major 
part of the space, where the interference is destructive and the wave function factually 
equals zero we then have to spend the computational resource only to check it. The 
dynamical diffusion swarm, in contrast, realizes the more general, Lebesgue integration 
scheme. In this scheme, the diffusion dynamics will concentrate diffusing particle in the 
areas of the constructive interference that gives the great saving of the computational 
resource. This is the fundamental advantage of the diffusion dynamics. We will see that 
the cost of this advantage is the non-uniform intensity of the diffusion on the element of 
space. This intensity will depend on the grain 5x of spatial resolution chosen in advance, 
in contrast with the usual diffusion, where there is not such dependence. 

We proceed with the definitions. We call a swarm a finite set S consisting of n point 
wise classical particles of the same type, each of which s & S possesses its coordinates 
and impulse p{s) G . The collective behavior method represents one quantum 

particle of the mass M and the charge Q by the swarm which samples s & S have 
the mass m — M/n and the charge q = Q/n. Members of this swarm we call samples 
of the particle. We assume that the total number of samples n is so huge that it can 
serve as the approximation of the continuous media. If we need to pass to higher and 
higher resolution there will be samples in each cube of the corresponding subdivision of 
the space. The dispersion of the speeds of samples will grow when the grain of spatial 
resolution decreases. It means that we will have the separate swarms for the different 
values of the grain. 

The methods of determination of particle depend on the specific problem, so that 
particles are not necessary elementary in the sense of theoretical physics. The definition 
of what to treat as a particle presumes the fixation of the typical length AX and time 
AT, so that if the size of a particle is much less than AX, then we treat it as a point 
wise, and the time period AT must not be less than the time of process we consider. We 
also assume that the typical medium speeds of movements are much less than some limit 
speed for the movement of real particles c. For example, we can treat an atom as point 
wise particle in processes with AX > 10~^m and AT > 10~^°s. If we make decrease 
the value of typical length and time, then for the right dynamical picture we need to 
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consider the other set of elementary particles, for example, we can consider separately a 
nucleus and electrons inside of atoms. Fixing AX and AT, we must determine smaller 

segments Sx, St, which represent elementary steps of the video film, though they are lesser 
than the typical lengths and times AX, AT for the more fundamental processes than we 
regard. The gap between these values can be about 10~^° that makes the distinguishing 
possible) . In the process with the fixed energy, the lengths and times depend on the mass 
of particles. The separation of particles by their masses in electrodynamics allows the 
consideration of only electrons, because the typical distances of nuclei movements will 
be 1800 and more times lesser. We then can regard the chosen values AX, AT as the 
relative size of an imaginary screen and the imaginary duration of the film, and 6x, 6t as 
the screen resolution and the time of showing of one cadre in our video film. We choose 
5x and 5t as large as possible that make our video film informative. After this choice we 
can determine, which particles are quantum and which are classical. For it we compare 
their typical actions a = M(AX)^/AT with Plank constant h. If a < h the particle 
should be treated as quantum, otherwise as classical. Wc will sec that in the method 
of collective behavior the passage from one type of consideration to the other means the 
change of swarm size, e.g., does not concern the introduction of the different types of 
dynamics. In view of the reserve we noted, and the choice of resolution we will be able in 
the preparation of the film decrease the values Sx and St in order to form the right image, 
for example, by the additional subdivision of this segments and substantiate the quality 
of the approximation to solutions of Shredinger equation. We assume that the space 
is divided to the equal cubes with the side Sx, and the time to equal periods St. 

We introduce some value of speed c, which we treat as the limit for the moving of 
samples. Segments of distance and time we choose such that Sx ^ cSt. It guarantees that 
at each step of evolution the values of magnitudes obtained by the averaging on cubes 
with the side Sx will vary slowly that is necessary for the asymptotic approximation. 

The density of swarm in the point x is given by the expression 

P(M) = ^. (13) 

where N{r, t) denotes the number of samples occurring in the instant t in the same cube 
with the point r. For the comparison with the solutions of Shredinger equation we would 
launch in these definitions Sx — > 0, which means that we consider not one swarm but the 
sequence of swarms with densities p„ and increasing n. We will not do this in the order 
to simphfy notations; instead we assume that it is always possible to spht additionally 
the division of space to cubes so that Sx will be lesser in the fixed frameworks. We will 
write p{x) — \^{x)\'^, which means that 

Pn{x)^\^ix)\\n^oo), (14) 

where the convergence will be uniform, without the special mentioning. This sequence 
of swarms, realizing the approximation to the density of the wave function solution of 
Shredinger equation, we call the admissible approximation of quantum evolution. 



27 



Our main aim is to define the behavior of samples, which gives the admissible approx- 
imation to the quantum evolution. 

The main requirement to the simulation of quantum dynamics via collective behavior 
is the following. 

• The swarm dynamics simulates the quantum dynamics so that in each instant t the 
quantum probability equals to the swarm density 

\^{x,t)\' = p{x,t) (15) 

in each point x of the configuration space. 

• Each sample has its own history, e.g., it preserves its individual identification number 
during all the simulation. The types of samples exactly correspond to the types of 
real physical particles. 

• The behavior of each sample is completely determined by its own state, the state of 
all samples in its vicinity and some source of random numbers. 

A swarm satisfying these conditions we call a quantum swarm for one quantum particle. 

We define the behavior of samples so that these conditions are satisfied. For this it 
would suffice to show that for any solution \l/(a;,t) of Shredinger equation we can move 
samples only locally, e.g., to small distances and thus ensure the equality f|T5|) in each 
time instant. 

The second rule means that we refuse from the complex numbers and at the same 
time want to include QED in our model in future. The last requires the locality of all 
interactions. The behavior of a sample is the rule determining the change of its state 
(impulse, momentum of impulse and the type) and the spatial position (spatial shift). As 
we know, the behavior cannot rest on classical mechanics. 

We define the quasi-classical behavior of samples called the dynamical diffusion mech- 
anism. The swarm with this mechanism will satisfy our conditions. 

We accept that each sample in each instant can either stand at place, or move with 
the speed c along on of the coordinate axes OX, OY, OZ. 

The reaction of exchange we call the following sequence of actions: a) the choice of 
the pair of samples a, (3, the distance between which is not greater than Ax, which speeds 
are mutually opposite: v{a) = —v{f3) and b) either the simultaneous replacement of their 
speeds to zero (if they are nonzero), or the acquirement by them the mutually opposite 
speeds of the module c directed along one of coordinate axes. 

The exchange does not change neither the summed impulse of the swarm, nor the 
summed momentum of impulse provided Ax is sufficiently small. We denote by N{r) and 
Ns{r) the set of samples in the cube with the point r and the set of samples with zero speed 
inside it, by N^{r), Ny{r), N^{r) we denote seta of samples from the cube with r moving 
along the corresponding axes in the positive direction, and by the analogous symbols with 
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the sign — in the opposite direction. By \g\ we denote the quantity of samples in a set g. 
We agree to denote the quantity of samples in a set by the same symbol, as this set, but 
with the replacement of by n. We call r - stationary each subset 5* C n(r), consisting 
of samples with nonzero speeds for which Ylaes'^^i^) ~ ^ ^^"-^ ^ maximal with this 

property. The number IS"! of elements of r - stationary set (which does not depend of its 
choice) we denote by s(r). Let c? > be the chosen constant so that the coefficient of 
diffusion is proportional to d, V{r)- is the scalar field proportional to the external potential 
energy with the constant coefficient of proportionality, grad V{r) = (V^(r), Vy{r), Vz{r)). 

We will consider only non relativistic swarms, e.g., such that ns{r)/n{r) is close to 
1 for all r. It means that the bulk of samples in each cube have the zero speed. This 
requirement is incompatible with the point wise approximation [T3] of the exact wave 
function by swarms for external potentials of the Coulomb type 1/r because the mean 
speed near the initial point goes to infinity. To have the asymptotic convergence [H] we 
must be able to choose the speed c as large as we need for the regular swarm of the number 
n. In the reality c cannot exceed the speed of light that establishes the natural limit to 
the accuracy of the swarm approximation of the solutions of Shredinger equation. 

We call the dynamical diffusion mechanism of the evolution the sequence of the fol- 
lowing actions with the swarm. 

• 1) The sequence of random exchanges with the uniform distribution leading to the 
distribution of speeds with the property s{r)/ns{r) = d for each point r. If n(r) is 
small, this equation must be true with the maximal accuracy. 

• 2) The ascription of the speeds to some samples from Ns{r), chosen randomly from 
the uniform distribution so that the signs of newly obtained speeds along each 
coordinate axes are the same and if f„(r) is the summed speed obtained by samples 
from the cube r along the axes u, u = x, y, z, then for such u the equation Vu{r)m = 
—Vu{r) is satisfied with the maximal accuracy. 

• 3) The change of coordinates r{a) of each sample accordingly to the law of uniform 
movement: r„ei«(a) = r{a) + v{a)At. 

• 4) The converting of V{r) accordingly to the new coordinates of particles. 

We do not specify the method of converting of the potential energy. We can do it by 
the Coulomb law or by the diffusion mechanism proposed in the work [?]. 

The swarm with the dynamical diffusion mechanism of evolution we call the dynamical 
diffusion swarm (DDS). This swarm is irreducible to the ensemble of point wise particles 
with the classical interaction. The point 1) says about two things: 

• there is the random force of repulsion or attraction, which preserves the summed 
impulse of the swarm (compare with [3]), and 
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• the averaging of speeds of samples goes with the accuracy determined by d (the less 
the value d is, the more exact averaging we have. 

For each time instant t if Ax is sufficiently small the density of swarm p{r, t) for any 
point r will not depend of the coordinate axes orientation. Really, let 5i be such that 

cSt -C SiX -C Sx, and v{r) denote the mean speed of samples in the point r, found by 
the averaging on the samples with coordinates ri : ||r — ri|| < 6ix. The total number of 
samples came in the unit of time from the vicinity of the point ri to the vicinity of the 
close point r2 will be proportional to the scalar product v{ri){r2 — ri)/\\r2 — Ti|p, which 
is independent from this orientation. 

We compare the defined dynamical diffusion swarm with the ensemble used in diffusion 
Monte Carlo method, which we call the stationary diffusion ensemble. The determining 
of the dynamical diffusion swarm state means the defining coordinates and speeds of its 
sample whereas only the density determines the state of stationary swarm because the 
mean speed of its sample is zero. Indeed at the step of evolution of stationary swarm 
for each its sample we choose with equal probabilities one of its shift: up, down, to right 
or to left, or it stays at place. This choice is independent from the initial speed of the 
sample and the resulting impulse is always zero. This model corresponds to the spreading 
of molecule of some substance in the media of the other substancej for example, of the 
color in the water. If we trace only the molecules of the color, the impulse will not be 
stationary: it passes to the molecules of the water. The stationary model of diffusion we 
can thus call the model with the friction. 

In contrast to the stationary swarm, the dynamical diffusion swarm evolves with the 
conservation of impulse in each interaction of pairs of samples in the exchange. Not only 
its density, but also the summed impulse in each cube of the space then determines the 
state of the dynamical swarm. E.g., the pair of functions p{f,t), p{f,t) then determines 
the state of DDS. 

However, the most important advantage of DDS in the comparison with the known 
ensemble methods is that each sample has its own history. It substantiates the name 
of sample. DDS allows the uniform consideration of a sample: from the classical and 
from the quantum viewpoints. To pass to the quantum consideration from the classical 
and vice versa we need only to replicate samples (or, conversely, to merge them). Here 
the swarm corresponding to the classical particle contains only one sample. We can also 
regulate the allocation of the computer memory in the simulation taking into account 
that the exactness of the quantum description grows with the total number of samples in 
the swarm. 

3.1 Differential equations for dynamical diffusion swarm 

To prove the appropriateness of the diffusion dynamics for the approximation of Shredinger 
equation we need to pass to the differential equations with the functions of real variables. 
Here the difficulty is that such equations will depend on the elementary length Sx. For 
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example, the intensity of the diffusion process is proportional to {6x)~^. It makes impos- 
sible to launch 6x to zero as in mathematical analysis, when we use it to the classical 
physics. The fixed grain of spatial resolution 6x determines the status of equations, which 
we will write. We must choose it such that the approximation of the density p = 
of the wave function by the density of the diffusion swarm within {6xY is satisfactory 
for the considered process. After the fixation of 6x we can consider the dynamical dif- 
fusion swarm of the corresponding intensity and the differential equations approximating 
its dynamics, which will turn to equivalent to Shredinger equation. 

The following pair of functions thus gives the state of the dynamical diffusion swarm 

p{t,r), p{t,f), (16) 

where p is the scalar function of the density of samples, p is the vector function equal 
to the summed impulse of samples in this point, which we define as lim P(r, dx)/{dx)^, 

dx >0 

where P is the sum of impulses of samples occurring in the cube around r with the side 
dx. Here we assume that dx can be done substantially lesser than the chosen value of the 
grain 6x, which defines coefficients of equations on p and p. 

The dependence of the equation on the grain 6x looks as follows. The summed impulse 
p{t, f) will vary slowly when f varies to the values greater than Sx. But its derivative ^ 
is huge: of the order l/{6xY, and varies quickly. E.g., the graph of function p(t,f) will 
be smooth if we observe it with the big grain Sx, but if we increase the resolution by 
decreasing the grain 6x, we see that it looks like a saw with sharp teeth. The sharpness 
of teeth will be the more, the more is the resolution l/6x, and is limited by the maximum 
speed of samples c only (compare with [3]). 

We express this requirement in the form: 

The intensity of the impulse exchange between the neighbor spatial areas 
must be much more than the intensity of sample exchange. 

This condition is important, we call it the non-relativistic approximation and write it 
as f <^ c. Factually, this requirement means the deviation from the classical physics in 
the interaction between samples of quantum particle, as the special character of impulse 
exchange. In the classical impulse exchange, we would obtain the oscillation equation 
with the second derivative on the time. The impulse exchange in the dynamical swarm 
has the peculiar character when samples standing at place simultaneously obtain mutually 
opposite impulses. 

In view of the isotropic property of the diffusion process the change of density p(r, t) 
in the time and its second derivative has the following expression through the integration 
on the surface S{r) of the sphere of radius 6x: 

a^pjr^t) - 3 f 9p{r,t) ( ^, I ^ (17) 

~9P~ - i;^^^ J -^^n[ri)as[ri). 

Sir) 
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These formulas follow immediately from the definition of density of samples and are true 
for any mechanism of the movement. 

We deduce the law ^a, of change of the summed impulse of the swarm in the small 
ball with the center in the point f, caused by the movement of samples along the normal 
vector a to the surface of this ball of the unit length. The deposit to the summed impulse 
of the small ball consists of three values: 

• Penetration of samples, obtained the speed in the exchanges through the small 
element of square. 

• Penetration of samples, obtained the speed from the action of the potential V. 

• Penetration of samples preserved their speed, e.g., by the inertia. 

It follows from the definition of the dynamical diffusion that these deposits are corre- 
spondingly, —I grad pd, —up grad Vd and gppd, where I,K,g is the intensity of the 
corresponding processes. By the appropriate choice of the unit system we can guarantee 
that g = 1. The following formula expresses the dependence of the coefficients on the grain 
of spatial resolution, which is necessary for the approximation of Shredinger equation: 

■ " (18) 



In view of the non-relativistic approximation we can reject the last summand, which 
is substantially lesser than the first two for the small 6x. It gives the formula: 

dp 

— ~ — / grad p — np grad V. (19) 
We thus obtain the equation on the density of the diffusion swarm: 

= - J I 9r(^d p-Kp grad V)n{r') dS{r'), (20) 

S{r) 

where coefficients /, k are given by formulas [T51 

We prove that the quantum swarm satisfies the equation [2D] that means the approxi- 
mation of quantum dynamics by the swarm evolution. 



3.2 Non- uniform intensity of the dynamical diffusion swarm 

The coefficient at the Laplace operator in the diffusion equation equals the intensity of 
the diffusion. The intensity determines the number of samples passing through the unit 
square in the unit time. For the simulation of quantum dynamics, we need the diffusion 
process with the non-uniform intensity. It means that the intensity must depend on the 
chosen grain 6x of spatial resolution. 
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At first, we take up the case when the potential is constant: grad = 0. The diffusion 
process with the non-uniform intensity can rest on the peculiar imaginary mechanism, 
which we call threads. We illustrate it on the following example. Let us suppose that all 
samples move not in all space but along some closed curve (thread) determined by the 
isomorphic injection of the circle to the space: 7 : — > R^. We assume that the process 
of exchange of speeds goes along this trajectory and samples always stay on it, with only 
change in their linear density. It is equivalent to the imposition of the holonomic bound to 
samples. We suppose that the linear density and the speed of samples are almost constant 
in each point of this trajectory. We consider one cube with the points from the trajectory. 
The stream of samples through its border does not depend on its side because there is 
only one trajectory. This form of the diffusion process has thus the intensity proportional 
because the quantity of samples penetrating in the unit of time into the cube 
with the side 5x, is independent of and the density is the ration of the quantity of 
samples to the volume. This example is not completely appropriate because many areas 
turns to be without samples at all, but it can serve as the approximate model when the 
total number of samples is small. 



3.3 Equivalence of quantum and diffusion swarms 

Here our aim is to prove that the sequence of diffusion swarms is the admissible approx- 
imation of quantum evolution. We have defined a quantum swarm as the satisfying the 
equation [TSl which evolution is reducible to local movement of samples. 

At first we prove that there exists a quantum swarm, e.g., that it is possible to satisfy 
[15] shifting samples locally. Then we prove that the mechanism of the movement of samples 
coincides with the diffusion, which gives the main result. 

We consider the quantum swarm, starting from Shredinger equation 

= -^Anr,t) + V,^,{r,mr,t), (21) 

which we can rewrite as 

for the real and imaginary parts \E'^, of the wave fnction ^. We take an interest only 
in the evolution of the density of quantum swarm, e.g. the function 

p(r,t) = (v^^(r,t))2 + (vl/^(r,t))2. 

We fix the value 6x and apply for the approximation of the second derivatives the 
difference scheme of the form 



9^\E'(x) ^ \l/(a; + Sx) + "^{x - Sx) - 2'${x) 
dx'^ {6xy 
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for each tine instant, treating the wave function as satisfying all the conditions for the 
accuracy of this approximation. Since the addition of any constant to the potential energy 
Vpot does not influence to the quantum swarm evolution, we can regard instead of Vpot 
the other potential of the form V = Vpot + ct) where a — —i^^^-i which eliminates 
the summand 2^(,t) in the difference schemes for the second derivatives on x,y,z (and 
thus gives the coefficient 3) after its substitution in Shredinger equation. To simplify the 
notations we introduce the coefficient 

h 1 

7 



2M {6x)^ 



Since we yet do not know the mechanism of movement of the samples in quantum swarm, 
we suppose that we either take some quantity of samples from each cube, or place them in 
cubes from some storage. We divide the quantum swarm evolution to so small segments 
of the duration St, that on each of them samples shift in the framework of neighbor cells 
only. If we prove that the diffusion mechanism gives the quantum swarm evolution on 
such a segment, it will be true for the whole evolution because our supposition about the 
exchange between only neighbor cells does not limit the generality. We also assume that 
cells differ one from the others only in shifts on Sx along the axes x, which also does not 
limit the generahty. We denote the centers of these cells by x and xi — x + Sx. Then due 
to our supposition about the exchange the summand \E'(x — Sx) in the difference scheme 
also disappears and on the considered small time period the following system of equations 
determines the quantum swarm evolution: 



fl(x) = -7^'(xi) + V{x)-^'{x), 



^j(x) = 7*'-(xi) - V{x)^'-{x), ^^^^ 

and the analogous system obtained by the replacement of x by Xi and vice versa. This 
system is applicable in the supposition that samples move from x to Xi. If they move 
from Xi to X the analogous system arises, with the replacement of x by Xi and vice versa. 
Without loss of generality we can choose the periods dt so small that on each of them 
movement goes only in one of these directions. In view of the complete symmetry, it would 
then suffice to consider the case of movement from x to xi, the second case is analogous. 
Shredinger equation then expresses the evolution process consisting of these two cases, 
where x and Xi can dispose by six different ways along three coordinate axes. Now by 
the period dt we mean the small period when the movement goes only from x to Xi. 
For this period we have 



^ = 2^*(x)(7^'-(xi) - V{x)^'Xx)) +2^''(x)(-7*^(xi) - V{x)'^%x)) 



27(^*(a;)^'-(xi) - W{x)'i/'{xi)) 



at ■ 



(24) 



It means that the loss of samples in one cell equals their return to the other that is the 
quantum swarm evolution satisfies the condition of locality. Now to compare its evolution 
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with the diffusion we find the second derivative: 



(-7^^(xi) + V{x)<i/'{x))<iJ'lxi) - ^'■(x)(7^''(x) - \/(a;i)^''(xi))] = 
272(^'-(xi))2 - 27l/(x)^"(x)^' (xi) - 2-f\^'{x)y+ 

27V(xi)^^(x)^^(xi) + 272(^*(xi))2 - 27l/(x)^^(x)^*(xi)- (25) 
272(^'-(x))2 + 27\/(xi)^''(x)^''(xi) = 

2^^{{w{xi)Y + C^K^i))^ - ((^'■(a;))' + (^*(^))'))+ 

27[(V^(xi) - + (^^(a;))^) + o(5x)], 

o(5x) = (^"(x)^"(a;i) + ^^(x)^*(a;i) - {{^'{x))^ + {^'{x)f)){V{xi) - V{x)). We now 
compare it with the expression for the second derivative of the diffusion swarm density 
found earlier, taking into account that in our case the exchange goes only between the 
neighbor cells along the axes x. Comparing with [20] and using [181 "we conclude that the 
second derivative of the quantum swarm density approximates the second derivative of 
the diffusion swarm density. 

If we choose as the initial state the state with the density of Gauss form, which 
coincides with the density of the ground state for oscillator, with the corresponding energy 
V = a(x^ + y'^ + z^) we will have dp/dt = in the starting moment in all points. By 
the proved, the second derivative of the diffusion swarm and quantum swarm densities 
approximately equal since the diffusion swarm serves as the good approximation of the 
quantum density on some interval AT. Including some potential slowly, we obtain the 
approximation for all quantum evolution in the limit of swarms when n goes to infinity. 

The swarm approximating one particle quantum dynamics depends on the choice of 
5x. After this choice we obtain the approximation of the wave function with the error of 
the order 5x provided 6t we choose as small as required. The intensity of the diffusion 
depends on the chosen grain, namely it is ^3^(^)3 • If we then want to decrease the step of 
spatial resolution, we must admit the more portion of moving samples in the unit of the 
volume. It means the uncertainty relation coordinate impulse: the dispersion of impulses 
of samples grows when 5x decreases. In any case, for the obtaining of the dynamical 
picture we must fix the grain of spatial resolution 5x. 

If the total number n of samples is limited we obtain the model of quantum dynamics 
with decoherence of the inbuilt type. Such a model we can generalize to the many particle 
case, so it serves the approximation of quantum dynamics in the standard Hilbert for- 
malism, (see ([?]). The appropriateness of this scheme for the numerical models follows 
from that it gives the direct probability space for Born rule, in contrast to Copenhagen 
approach, which takes Born rule as an axiom. 
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3.4 Restoration of wave function from the dynamical diffusion 



We have solved the problem of approximation of the dynamics of quantum particle density 
by the special diffusion process with the non-uniform intensity. A pair [16] determines a 
state of the dynamical diffusion swarm. Such a pair does not use the notion of complex 
numbers, which generate quantum interference in the standard formalism, and it does 
not give the beautiful differential equations of Shredinger type for p and p. Moreover, 
the mechanism of the dynamical diffusion we introduced for the simulation of quantum 
evolution differs from the classical processes (as oscillations or heat transfer) in that its 
intensity depends on the chosen grain of the spatial resolution. We put up with it for the 
sake of the main: the economy of the computational resources required for the simulation 
of quantum dynamics. 

Now for the completion of the picture we have to solve the converse task: to chow how 
to restore the conventional wave function ^ given a state of the diffusion swarm [161 For 
this we address to the equality [211 cind substitute in it the expression fro the wave function 
through the density: \E'(r) = -\/ p{r) exp(i0(r)). We need to find the phase 0(r) of the 
wave function. We note that since only relative phase between points has the physical 
sense, we can fix some point r and consider a phase in some other point ri relatively to 
r. If ri close to r, the equation [Ml gives 



where the path 7 goes from r to ri. This definition explicitly depends on the choice of 
the path 7, hence we have to prove its correctness, e.g. the independence of this choice. 
Since the phase is determined within an integer multiplier of 27r, the different choices of 
the path can result at most in the change of the phase to this number, as for the excited 
electron state in a hydrogen atom with the nonzero momentum (for example, 3(i). We 
show that the integration of the speed v of the swarm along the closed contour preserves 
its value in the time the more exact the less is the grain of spatial resolution 5x. If then 
the definition (1261) was correct in the initial time instant, it preserves its correctness for 
the further evolution as well. 

We thus consider the derivative of the integral of the speed of swarm along the closed 
contour 7c. Applying the formula f[20l) in view of dp/dt = p dv/dt, we obtain 



swarm 




that leads to the following formula for the determining of the relative phase: 




(26) 



7 




(27) 



36 



The first summand gives zero after tlie integration along tlie closed contour because it is 
grad In p, the second summand gives zero by the analogous reason. 

It would now suffice to check the correctness of the definition (l26l) in the starting 
moment, which we can do for each concrete problem. In the case when the wave function 
of the initial state comes from the ground state in Coulomb field where w = 0, the 
correctness follows form the proved because here there is no phase shift to 27rfc. If to find 
the initial state in the considered problem we need to start from some excited state with 
some storage of the phase, we need to check the correctness of the definition for this state 
first. 

3.5 The case of many quantum particles 

We consider swarm representations of n quantum particles 1, 2, . . . , n, where Si, S2, ■ ■ ■ , Sn 
are swarms of samples corresponding to their states l^&i), 1^2), • • • , l^n)- If we consider 
the ensemble consisting of all these samples, it represents the non-entangled state of the 
form |\E'i) (3) 1^2) 'S) ■ ■ ■'S) l^n)- However, to represent entangled states of the form 

Yl ^iiJ2,...jJjl,i2,...,in) (28) 

We need to introduce the new essential element in the method of collective behavior. It 
is bonds between samples from the different swarms. Any basic state ji we can treat as 
the coordinates of some particle i in the corresponding configuration space. The repre- 
sentation of the wave function in the form [28] means the existence of bonds connecting 
points ji,j2, • • • , jn in the same cortege. The relative quantity of such type of bonds (their 
number divided to the total number of bonds) is |Aj^j2^...j^ p. 

We assume that bonds connect not spatial points but samples of the different real 
particles. These bonds we write as corteges 

S = {Si,S2, . . . ,Sn) (29) 

where for each j = 1,2, . . . ,n sj G Sj. The wave function |$) then acquires the form of 
the set S of corteges s so that for all j = 1,2, . . . , Sj E Sj there exists exactly one cortege 
of the form [29l Each cortege plays the role of the so-called world in Everett many world 
interpretation of quantum theory!^ We treat this cortege [22] as the probe representation 
of the system of n particles and assume that all interactions go inside one cortege only, 
whereas the state of the real system resulted from the interference of amplitudes of all 
corteges occurring in this spatial cell, or resulted from the other process replacing the 
interference of amplitudes. In the case of one quantum particle, we saw that the exchange 
of impulses in the diffusion dynamical swarm could serve as such a process. We now 

Strictly speaking, for this we must treat the times t as the part of corteges (see [4]). In the non- 
relativistic theory there is the single time t. 
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study the generalization of this process to the case of n particles. We call 5" the swarm 
representation of the system of n particles. 

The density of swarm S is given by the formula 

p,(r„r2,...,rO= ,lim ^-;r^:;;;;' , (30) 



dx — »oo (dx] 



where iVri,r2,...,r„, dx is the total number of corteges locating in 3n dimensional cube with 
the side dx and the center ri, r2, . . . , r„. 

If the wave function |$) is the tensor product of one particle wave functions: 



= (g) m 



1=1 



then the corresponding bonds we can obtain by the random choice of samples from the 
uniform probability distribution sj G Sj for each j = 1,2, ... ,n, which thus form each 
cortege si,S2, . . . , s„. With such a choice of corteges we obtain that the density of the 
corresponding swarm satisfies Born rule, which we can read for swarms in the form 

E\(^m\' = ^ (31) 



where D C -R^", Nf,^g is the total number of corteges located in the area D. For the 
entangled state |$) this choice of corteges for the set of swarms S does not give us 
the condition [3TJ We thus have to take [31] as the definition of the corteges in S. To 
determine the swarm we must also define the speeds of all samples, namely, to generalize 
our construction to the case of n real quantum particles. 

Let ^'(ri, r2, . . . , r„) be the wave function of the system of n particles, 
\1/ = |\l/|ea;p(z0(ri, r2, . . . , r„)) is its Euler representation. We denote by gradj0(ri, r2, . . 
the gradient of taken on coordinates of the particle j, where j G {1,2, . . . ,n} is the 
fixed number. For n particles the formulas of the relation between the wave and swarm 
representation of the quantum state have the form 



|^(f)| = 

0(r) = / k{dxYv ■ d-f, (^32) 



7: ro- 



V = a{dx) grad 0(r), 

where f denotes ri, r2, . . . , r„, grad denotes (grad^^ , grad2 , . . . , grad„), and 7 is the path in 
3n dimensional space. The rules [32] are sufficient for the determining of the swarm given 
a wave function, provided we agree to join samples into corteges independently of their 
speeds. The microscopic mechanism of the swarm dynamics then acquires the following 
form. The exchange of impulses between two corteges of samples: s = (si, S2, . . . , Sj, . . . , s„) 
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and s — {si, S2, . . . , Sj, . . . , Sn) is the impulse exchange between two samples Sj and Sj 
provided s and s belong to the same spatial cube in the configuration space i?^" for n 
particles. Here j is chosen randomly from the uniform distribution. With this definition 
the reasoning we earlier carried out for one particle can be repeated literally, and we 
obtain that this microscopic mechanism of impulse exchange for n particles ensures the 
approximation of n particle quantum dynamics within the accuracy of the order dx^'" for 
the definition of the wave function on the fixed time period. 

The method of collective behavior gives us the basement for the economical modeling 
of the real quantum evolution, including the unitary evolution and decoherence. 
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4 Lecture. Heuristic of collective behavior 



We have established that the simulation of quantum evolution by swarms with samples 
joint in non-overlapping corteges gives the approximation of the unitary quantum evolu- 
tion provided the quantity of samples of each particle goes to infinity. Factually, the model 
of a quantum system with many particles we represent in the form of two-dimensional net 
of the form shown at the picture. Here in each column stand samples of one quantum par- 
ticle and in each row stands the quantum world, to which from each real particle belongs 
exactly one sample. Quantum exchange of impulses goes on each real particle between 
close worlds, and the potential energy acts in the framework of one world only, e.g., each 
sample is subject of influence from only those samples, which belong to the same world 
with it. 

The dynamics determined by this net differs from the unitary only to the extent of 
the total number of corteges (equal to the quantity of samples of each particle) is not 
infinite. If the quantity of particles grows the rapid divergence of the net dynamics from 
the unitary, determined by Shredinger equation takes place. It is just the absolute model 
of decoherence, when the lack of classical memory for the keeping of quantum states 
plays the role of decoherence. However, it is not sufficient for the building of modeling 
algorithms. We must somehow determine how many samples of real particles do we need, 
e.g., what is the real value of n. We can treat the determining of the real value of n in 
the nature represents the purely physical problem admitting the exact solution. I do not 
know it. Though we can assert that the value of n can come only from the observing of the 
dynamical scenarios built by the net. Moreover, this value can depend on the simulated 
process in the following sense. Processes with the small length of corteges Si, S2, . . . , 
can allow the bigger values of n, whereas the process with the large m conversely presumes 
the small values of n. In any case, the capability of the dynamical scenario influences to 
it much more than the known physical constants. 

What new the swarm representation gives for the description of many particles compar- 
atively to the standard quantum formalism? The new possibility is that the entanglement 
obtains the direct representation in the form of corteges Si,S2, • • • ,-5^, which consist of 
separate samples of all particles connected by the bonds. We call the bonds the objects, 
which connect the samples to the corteges quantum worlds. Bonds have no analogs in 
the standard formalism, which also has no samples. The standard treatment defines the 
entanglement as the impossibility to represent the state vector as the tensor product of 
one-particle states. It is no way to work with this definition from the contrary but the 
matrix algebra, where decoherence is the foreign element, which thus does not admit any 
analysis from the first principles. Nevertheless, it is reasonable to treat the entanglement 
as the physical resource, as particles themselves. It means that we must ascribe to the 
entanglement some formal object making possible to manipulate with it directly. Only 
constructivism gives this possibility. In the heuristic of collective behavior, we can in- 
troduce rules of the break of a long bond when the (conditional) burden to this bond 
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Figure 1: Net representing the system of m quantum paritcles 
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exceeds some threshold. The other form of the work with bonds is their selection. The 
selection of quantum states has the following form. Corteges worlds we can join to the 
separate groups Fi, . . . , F; so that any group will contain corteges close to each other 
in the configuration space for m particles. We then can cancel groups with small quantity 
of samples and redistribute bonds between the rest samples accordingly to the closeness 
to worlds in large groups. This principle of quantum selection has the interpretation in 
terms of constructiveness of the amplitude interference. 

We consider d samples each of which carries the amplitude of the module a, but with 
the different phases. If these samples stay near each other, their phases are close and the 
deposits of these samples in the resulting probability are approximately (da)^, because we 
must add just their amplitudes. If these samples are far, their deposit to the probability 
is d{aY, which is in d times lesser, because now we have to add probabilities (phase are 
distributed randomly). 

This reasoning is valid for the electromagnetic field as well. Given an ensemble of 
atoms, which can emit photons, we have two possibilities: 

• The change of state of the atom, emitting the photon makes possible to identify it 
among the others atoms of this group. 

• The change of state of the atom, emitting the photon redistributes among all atoms 
of the group that makes impossible to identify this atom. 

Let ip be the amplitude, which the emitted photon carries, ^ be the state of the group 
of atoms. We ascribe the emitting atom the index j and assume that amplitudes Xj of the 
emission of each atom j arc approximately equal. Let the zero vahic j = corresponds to 
the case when the photon is not emitted. In the first case, the state of the whole system 
of atoms and photons is determined by the common wave function of the form 

E^^-|^^-)l*^)' (33) 
j 

in the second - 

J]A,#,)|*) = |V',ene.a;)|*, (34) 

3 

where \ipgenerai) = S'^jlV'j) is somc photon state. We estimate the probability of the 

3 

photon emission in these cases. In the first it is 

3^0 

in the second 

lE^^f- (36) 

3¥=0 
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Since all Xj for j ^ are approximately equal, the probability of the photon emission in 
the first case is much lesser than in the second: the coefficient is the total number of atoms. 
This qualitative reasoning belongs to Feynman (see [2]). Its construction interpretation 
rests on the notion of bonds between samples of atoms and photons. In terms of collective 
behavior, the redistribution of the return impulse between atoms in the photon emission 
goes through the bonds between samples of atoms. In standard terms, atoms must be 
in the entangled state; however, the detailed description of the redistribution is hardly 
possible, because just such problem requires a quantum computer. The constructive 
heuristic allows the simple reasoning. The entanglement means the presence of stable 
bonds between the samples of separate atoms. The redistribution comes from the kind 
of elastic forces generated by bonds in this world. This property of elasticity is not 
the additional element of the reality, but the main property opf bonds; for example, it 
ensures the interference properties of a molecule which wave function has GHZ type: 
^ Aj|\E'jj) |\E'j2) . . . I^E^jj.), where |\E'j) are states of its atoms. 
] 

The existence of bonds between atoms thus substantially increases the probability 
of the photon emission. We study the difference between the separate photons and the 
electromagnetic field. This difference is that the field is the ensemble of photon samples 
connected by bonds. In the standard formalism the photon wave function in our example 
{ipgenerai) represents (after the appropriate renormalization) the intensity of the classical 
electromagnetic field. In the first case, when we have the separate photons, we can speak 
about the field only conditionally, e.g., after the detecting of the concrete photon and 
this filed if very small. In the second case, when the return impulse from the photon is 
redistributed, we have the solid value photon amplitude, which reveals classically. This 
is the standard description of the photon radiation. We see its duality: in the case of 
the field we cannot speak about the probability, but only about the field, measurable 
classically, whereas in the case of separate photons we speak just abut the probabilities. 
The standard approach makes impossible to use the same language for the description of 
the emission in the cases of separate photon and the field. This is the feature of standard 
formalism. 

In the constructive formalism the both situations with the separate photons and the 
field have the uniform description. The field differs from separate photons only in the 
presence of bonds between samples. What we earlier called the elastic properties of these 
bonds guarantees the possibility to observe the field by classical means. It is important 
to stress that bonds are not material objects; it is only the mathematical abstraction 
belonging to the constructive version of many particle quantum mechanics. It is thus 
senseless to ask questions about their mass, etc. Analogously, what we call the elasticity 
of bonds is not the usual elasticity obeying Guk law. It is the heuristic trick only, which 
helps to solve the main problem of constructive quantum theory: the creating of the 
working model of dynamics. 

This situation shows why the detection of the entanglement of quantum states of 
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a complex system can be difficult. The entanglement is not the characteristic firmly 
connected with the fixed particles, but can pass to the others particles. Let us consider 
the system of three qubits in the standard formalism, when they are in the states ;^(|0) + 
|1)), |0), |0) correspondingly. The sequential operations CNOT at first on the first two 
particle, then on the first and third give the state of the form -^(|000) + |111)), in which 
all three qubits are entangled. If we now apply the next operator CNOT to the first and 
second qubits it results in the entanglement of the first and the third qubits whereas the 
second will be non-entangled with them. The entanglement can thus easily pass to the 
other particles. In the various interactions in the complex system, this situation is the 
rule and it leads to the difficulty in the detection of the entanglement. In any case, in the 
complex systems the methods of detection of the entanglement working for ions in the 
trap or for biphotons do not work. 

Nevertheless, the general scheme of the detecting of entanglement is also apphcable 
to complex systems, provided we have many identical copies of their states. This scheme 
is the quantum tomography. Given a state of two qubits Si and ^2, we must determine 
is it entangled or not. For this we need to prepare this state twice and measure it in 
the basis E — and in the basis E' — E[<^E'2. If Ei coincides with E2 and 

E[ coincides with E2, the results of measurements will be correlated. For example, we 
can thus distinguish EPR state ;^(|00) + |11)) from any mixed state with non-entangled 
components. We can vary this scheme and fit it to more complicated systems ^i, 5*2, we 
also can choose two basic systems Ei and E[ differently. The main is that in the case 
of the same choice of these basic vector systems the results of measurements must be 
correlated. In the constructive quantum mechanics, this property is the characteristic of 
bonds between samples of the different particles: we must take into account the reaction 
of the whole quantum world in the interactions between samples. Of course, it does not 
mean that the reaction of samples connected by bonds must be equal; it may be opposite. 
The main is the presence of some correlation in the reactions of all samples, belonging to 
this quantum world. 

The various examples of such correlations are well known. We consider two living 
things and their reaction to the same external influence. If this influences are graduated 
so that we can compare the reactions of the things to these influence, these reactions will 
be the more similar, the closer these things are. We can introduce the degree of closeness 
of hving things by their DNA closeness. The degree of closeness of their reactions will 
be then in the direct dependence of the degree of their affinity. This closeness may be 
almost absolute, as for the one-ovum twins. It brings the analogy with the entangled 
quantum states in which different objects may be. However, the quantum entanglement 
is not reducible to classical correlations. It follows from the violation of the so-called 
Bell inequalities, which we discuss later. The peculiar sense of quantum correlation is 
that the nature behaves as the constructive model, in which the user cannot interfere 
into the created video fllm, but must observe it from the beginning to the end. Its free 
will permits only to order to the system to create such a video film. In the process of 
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building of it, the communications between spased points takes place, which the user 
can treat as the instantaneous communications. It is necessary to show the quantum 
correlations in the film. These communications belong to the administrative part of the 
model, inaccessible for the user. In particular, the user cannot apply these communications 
for the transmitting the information he created to the far distance, so that the principle 
of relativism remains valid. If we assume this limited treatment of the free will as the 
natural, no obstacles for the constructive representation of the entanglement in the form 
of bonds remain. This representation is the direct analogy of the quantum entanglement 
with the classical connections; we must accept this treatment if we want to enter the 
world of complex systems. 

In the standard formalism nothing stays behind this analogy and it is the superficial 
similarity. Indeed, in any reasonable specification of the term hidden parameters it is 
possible to prove their absence in quantum theory (one version of it is Cohen Spector 
theorem, there are the others). It means the following. In the standard formalism there is 
nothing but the wave function that affects the result of the measurement. This postulate 
of quantum mechanics is called the theorem about the absence of the hidden parameters. 

The theorem about the absence of hidden parameters has the sense only in the standard 
formalism. Namely, the sense of this theorem is that the space of elementary events for 
quantum probability is inaccessible for us. It the constructive quantum theory we always 
can use explicitly elementary events operating with corteges of samples, which compose 
quantum worlds. It gives us not only the direct form of elementary events for quantum 
probability but also the possibility of the immediate manipulation with non-local states 
in complex systems, which role is yet to be elucidated. 

4.1 About fermionic anti symmetry of wave function 

We consider the representation of the fermionic symmetry of the wave function in terms of 
collective behavior. Let denote one particle wave functions, and |i) denote its spatial 
positions. Accordingly with the requirement of the anti symmetry of the wave function of 
the system of the independent identical fermions it must have the form |^) = det 
where in the expansion of the determinant all multiplications we treat as tensor products. 
Such states |^) we take as the non entangled states of the system of identical fermions, 
whereas all other states can be expanded to the rows on these functions taken for the 
different choices of one particle From the view point of Hilbert formalism, such 

states are entangled. However, in the representation of collective behavior we can easily 
make them non-entangled. Really, imagine that the supports (area where the function 
is non-zero) of them do not overlap. Then the determinant becomes the product of one- 
particle functions. We can make the supports non-overlapping by the decreasing of the 
grain of spatial resolution dx. In the swarm representation, two swarms cannot coincide 
even in one point due to the Pauli principle, and the determinant we can then replace by 
the ensemble of two swarms. By the way, it makes possible to specify the particles, which 
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are identical from the canonic viewpoint. This specification follows from the completely 
different supports of their wave functions. I stress that in the collective behavior the 
fermionic character of particles belonging to the same type follows from Pauli principle 
and has the representation as the non-overlapping supports of their wave functions (not 
intersecting swarms). 

4.2 Constructive viewpoint to quantum electrodynamics 

We look how to make constructive the formalism of quantum electrodynamics. The 
ordinary (non-relativistic) quantum theory deals with the wave function dynamics 
where the physical time t is the common parameter; the case of relativism is the different. 
Here we have a state |^'), which is in the following relations with the time. There is 
the so-called internal time of quantum system t, and there is the external time r. The 
evolution accessible for the user observation is the sequence of the unitary operators in 
the space of states 7i of the form 

|^,=o) |^r=l) . . . ^ 1^.=^) (37) 

called scatterings. This is the chain of the quantum state evolution in the external time 
r, which periods we denote by integers. We choose them so that this sequence gives us 
the realistic picture, for example, of the chemical reactions. Here all states \'^T=j) = |^) 
have the form 

m = J2>^^\^^M (38) 

i 

e.g., is the linear superposition of basic states \xi)\ti). Each of such states means that the 
considered system is in the point Xi = {x},xj, . . . ,x") of the configuration space in the 
moment ti = {tl,tj, . . . ,t") of the internal time. Here xj is the coordinate of the sample 
j of the real quantum particle i in our system, is the time instant corresponding to 
this sample. The internal time of the system in QED is not thus the unit parameter, 
but something similar to the spatial coordinate: each sample has its own time. The 
presence in the expansion just the time cortege tj follows from the conventional rules of 
the operations with tensor products of states |^'') and |^^") for the different particles, 
when their state in a non entangled case has the form ^ \ '^'), and in an entangled 
state is the superposition of these states. Non relativistic states, in which the time is 
common is factorized as the tensor multiplier, can be only spatially entangled, whereas 
the relativistic states have no common internal time, and states can be entangled on the 
internal time as well. 

It is possible to prove (see [4J) that the law of energy impulse conservation in the 
fundamental interaction between charged particles and photons gives the presence of the 
common internal time for the interacting particles within the chosen grain. It means 
that inside each quantum world for close samples the time must be approximately the 
same because they must interact. Analogously, for the samples of the closed quantum 
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Figure 2: Network representing scattering 
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worlds (here the closeness include the times as the coordinate in the coordinate in the 
configuration space) exchanging by impulses the difference in the time must be in the 
framework of the grain. Wc obtain the following picture for the relativistic net of QED 
corresponding to one state Here the representation will be not two dimensional, as in 
the non-relativistic case, but three dimensional, because the internal time t now belongs 
to the configuration space and it has the same detailed structure as the space coordinate 

X. 

What means now the selection of quantum states in the evolution in the external time 
r? On each passage from r to r + 1 we select states of the form \xi) \ti). In is not the selec- 
tion of ordinary wave functions distributed in the space. This is the selection of dynamical 
scenarios, e.g., wave functions distributed in the space-time. How this selection goes? Its 
simplest form is the grouping of states by their closeness in the configuration space-time 
and the selection of groups with the maximal density. It is the simple generalization of 
the corresponding principle for the non-relativistic case. In this selection, factually bonds 
of two types are in use. 

• Bonds, joining samples of the different particles, e.g., bonds of the type we already 
discussed. 

• Bonds, joining one sample taken in the different instant of its individual time. 

The second type of bonds means the presence of individual histories of each sample. 
How these bonds on the internal time work in the modeling? They have the certain role. 
We select not only close states, but close dynamical scenarios. Such scenarios may be the 
attribute of concrete samples of real particles. Their regrouping in the selection in the 
genetic algorithm means the recombination of these attributes. The both these types of 
bonds belongs to the planes tti or 712 at the picture. 

Are there more complex methods of selections than the grouping of scenarios by their 
closeness? In the other words, is it possible to introduce bonds corresponding to sum- 
mands in the wave function of the form |a;2) [^2) for the different times ti ^ t^^. 
Such bonds would have the form of slanting lines in the picture, e.g., non-parallel to the 
planes tti or Ti2- The introduction of such bonds may come from the new methods of the 
dynamical scenarios selection; hence, this question remains without answer. 

The constructive consideration quantum system with many particles in QED leads 
us to the work with dynamical scenarios, each with its internal time t. It includes the 
selection of such scenarios. The surviving of certain scenarios replaces the surviving of 
the certain states in non-relativistic models. In these scenarios participate many samples 
of real particles, which can be in the scale of external time for a long time r, which gives 
the observed phenomena of the microscopic character. 
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4.3 How can we create quantum computer and why do we need 
it: realization of collective behavior on computational net- 
works 

The method of collective behavior gives the principal way to overcome the problem of 
quantum computer, e.g., the exponential growth of computational resources needed for 
the exact modeling of the unitary dynamics in the full Hilbert space of states. Here we 
treat decohcrcnce as the fundamental factor, represented as the limitation of the classical 
memory. It makes possible to simulate biochemistry on classical computers. 

However, in the heuristic of collective behavior we meet the difficulty as well. It 
is the necessity to transform the information along bonds joining samples of the same 
quantum world. This difficulty may seem empty to the pure mathematicians, but this is 
the real problem. The scientific logic joins samples not in quantum worlds but in swarms 
corresponding to real particles. It makes unavoidable the situation when the different 
groups of programmers work at the different groups of real particles. If we join real 
particles in some non overlapping ensembles Fi, . . . , F^, the evolution model of each of 
them Fj will work at the corresponding processor Pj. Bonds joining samples inside each 
quantum world then become the bonds between these processors. 

What kind of information can these bonds transmit? At first, it may the simple syn- 
chronization as the type of classical correlation. This kind of information is similar to the 
open key distribution in quantum cryptography. Let the different processors Pi and P2 
have to generate the common code, which will then work at them. This problem arises 
in the modeling of the specific for the living thing proteins gamma globulins allowing 
the identification of own tissues by antibodies. For this kind of problems, we can use the 
methods of one particle cryptography, including its quantum version. However, if the syn- 
chronization is included in the complex process of local computations on both processors, 
the preferable choice for it will be bi-photonic mechanism. This mechanism causes no de- 
lay in the computation. The idle time of processors coming from the passage of the signal 
from Pi to P2 may be of the order of hundreds microseconds if these processors are in 
the different towns, and it can valuably brakes the work of the whole net. This argument 
touches the organization of computations only, not the modeled process themselves; I do 
not assert that in the real systems the synchronization is carried out by this mechanism. 

The other example illustrating the importance of synchronization is the process of 
penetration of a molecule through the cellular pore. The cellular membrane has pores, 
which can be either in open or in closed state. Penetrating molecules transport important 
substances, and at the same time, the membrane must be impenetrable for the other, 
harmful substances. The pore thus must open just in the moment when this transport 
is ready for the penetration. The pore itself is the complex mechanism, as the protein 
molecule, it is thus non-reasonable to charge the same processor by the simulation of 
the internal evolutions of these two systems. The simulation of the whole process thus 
requires the synchronization of evolutions for the obtaining of the desired effect. 
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Figure 3: Penetration of transport molecule through pore in cellular membrane 
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I speak here about the instantaneous classical correlation of the different processors 
by pairs of photons. The usage of entangled photon pairs promises the more advantages. 
If the forming of hydrogen bonds between two complex molecules belongs to the modeled 
process, for the simulation of such quantum particle belonging to each system it would be 
useful to have the real quantum model of this bond in both parts of the computational 
process. 

Of course, we can simulate each quantum correlation by classical means if admit the 
possibility of the signal to travel form one processor to the other. The drawback of 
this ideology is that it contradicts to the basic requirements of the computational net 
architecture. These requirements set the general form of control shown at the picture. 

The central processor sends signals to the different local processors, each of which 
models the corresponding subsystem in the whole system. For example, CP decides that 
it is the time for the synthesis of some polymer A of the specific activity in some subsystem 
and simultaneously the other polymer B, suppressing (or intensifying) this activity in 
the other subsystem. CP sends the package of signals to the both subsystems, then 
immediately switches over to the next task (for example, the simultaneous synthesis of 
the other pair A' and B') . What happens if subsystems themselves begin to send signals 
to each other? Let we be given m different subsystems, ach of which controlled by its 
own processor. For the right addressing of mutual signals between local processors (about 

ways to choose them) we need to charge CP with the addressing of these signals. CP 
then will have to wait the time cD, where D is the distance between local processes, c is 
the speed of hght before the switch over to the next task. If there are many tasks (in the 
real model of biochemical metabolism there are very many such tasks) this delay results 
in the fatal braking of the model. This is why the idea of the direct exchange of signals 
between the different local processors is not good. We thus come to the necessity to use 
one side control from CP, which realizes it by sending of signal packages to the different 
pairs of processors. CP also receives signals from the local subsystems, but not directly, 
by the chain of connected processors. This form of the organization of a living system 
model is more appropriate because we maximally relieve CP from the routine tasks of the 
metabolism maintenance. 

The architecture of the computational net must then have the form shown at the 
picture, where the general control realizes CP. It communicates with the bulk of local 
processors so that the signal form the remote subsystem comes to CP through the chain 
of intermediate subsystems only. 

The collective behavior heuristic represents the entanglement of quantum states by 
the special object bonds between samples of the different real particles. We suppose that 
just manipulations with bonds makes possible to model complex systems, because this is 
the only way to work with entangled states. One processor will process bonds between 
samples belonging to the same subsystem. Bonds between samples from the different 
subsystems require the special attention. These bonds we call remote. 1 stress once again 
that we speak not about the real distant subsystems but about the distant local processes 
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Figure 4: Structure of complex system model 
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in the net. There are seasons to treat the remote bonds as important for the obtaining of 
the right scenarios of the complex system dynamics. Hence we must charge CP with t he 
work with the remote bonds. It can act purely classically, and send signals to the chosen 
local processes. 

There exist situations where CP using entangled bi-photons has the advantage over 
CP with classical signals. We consider the following abstract problem. Let we have to 
synthesize two molecules of polymer 

Ai, A2, . . . , Am 
Bi, i?2, . . • , Bm 

consisting of mono blocks of the two types: a or b (see the picture). Each mono block 
has the lateral surface, which can glue if it adjoins to the similar surface of the other 
mono block, and the internal surface, which can glue to the other mono block only by the 
special boll in its center (see the picture). This construction ensures the non-symmetric 
gluing when two mono blocks are shifted cither in the same direction, or in the opposite 
directions. Two glued mono blocks we treat as firmly fastened. Mono blocks are linked by 
the intermediate linear segments, which can stretch or shrink. Here each mono block can 
be shifted along the chain to some dx either to right or to left. The synthesis goes by the 
sequential addition to the existing chain the new mono block, which comes to the point 
of assembly first (mono blocks in the storage are in the chaotic movements). It is possible 
to shift a mono block either ahead, or back along the chain to the small distance dx. We 
denote the shift ahead by +, the shift back by — . Each segment Aj or Bj is the pair of 
the form c, s, where c e {a, b}, s G {+, — }. Let the assembly of the polymer A goes in 
one point, and the assembly of the polymer B - in the other point so that their models are 
located far one from the other (for example, they are located in the different countries). 
The problem is to organize their synthesis so that the quantity of critical impositions in 
the matching of the chains A and B is as small as possible. We call the critical imposition 
such imposition, in which the mono blocks are not firmly fastened. 

Such problem could appear in the modeling of the simultaneous synthesis of gen and 
its anti-gen in the different cells. Of course, we could organize the simulation by the 
information exchange between the two points of assembly at each step. However, if the 
corresponding processors are far one form the other this exchange would substantially 
impedance the modehng because we have to wait when the signal travels between com- 
puters at each step of the process that reduces the efficiency of the using of fast com- 
puters. It is possible to organize the synchronization by the classical correlated photon 
pairs. However, the using of the entangled bi-photons in this problem gives the valuable 
advantage over the classical schemes. 

Since for the minimization of the total number of critical impositions we intend to use 
the state of EPR type |^) = -i^dOO) -I- 1 11)), it is necessary to substantiate that it gives 
the advantage over the classically correlated photon pairs. For this we recall Bell inequal- 
ity, valid for all classical correlations. Let we be given four classical random variables 
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b2,bi,ai,a2, taking values from the set {1,-1}, for their mean values Bell inequality 
takes place: M{\aib2 + 6162 + 0102 — &ia2|) < 2, which follow from the representation 
0162 + bib2 + 0102 — 6102 = ai(&2 + 0-2) + bi{b2 — 02) because one of the brackets is zero 
and the module of the other is 2. Wc accept the agreement. The lower index denotes the 
number of the point of assembly: 1 or 2. The letter a or 6 denotes the type of mono block, 
and the value of the corresponding value equals +1 if we shift the mono block ahead, and 
— 1, if it shifts to right in the assembhng. The result of the assembhng in the both points 
is determined provided for each lower index 1 or 2 we have at first, the letter a or b, and 
at second, the sign of the shift + or — . The letter a or 6 is always determined by the type 
of the mono block, which at this moment randomly occurs in the closest distance from 
the assembly point. We thus can control the assembling only choosing signs of the shifts: 
+ or — in the both assembly points. To guarantee that the communication between them 
docs not decelerate the modeling of the assembling process, we need to choose two sings 
in both points simultaneously. In the case of classical correlation of these choices, we 
would have Bell inequality. For each step of the assembling process we introduce the 
index of criticality as Cr = +1 if it gives non critical imposition, and Cr = —1, in the 
opposite case. We are interested in the total number of non-critical impositions in the 
whole polymer NonCr. Our aim is to make it maximal. For one pair of mono blocks 
we have NonCr = |(1 + Cr). Because each combination aa,ab,ba,bb occurs with the 
probability 1/4, for the mean value M{Cr) of the index of criticality we have 



which by Bell inequality gives us for the mean value of the number of non-critical impo- 
sitions the inequality M{NonCr) < |(1 + |) = | = 0.75. 

In the case of quantum EPR state of bi-photons the situation is the different. There 
is no the common probability space for two point of assembly, where a and b with indexes 
are random variables. Of course, there exists the common probability space of elementary 
events P for quantum probability, but in this space we cannot treat a and b as random 
variables, because this general space P consists of non-local events. It means that here 
we have not Bell inequality and must compute the probabilities directly, using Born rule. 
Let in each assembly point stands the corresponding photo detector, which can obtain 
orientations corresponding the desired observables a and b. For the first and the second 
points, these observables have the form 



correspondingly. We agree that the type of the current mono block in the assembling 
determines the corresponding position of the detector in each point, and the shift sing 
of the mono block equals the value of the corresponding observable. Again, since all 
combinations of the monomer types aa, ab, ba, bb have equal probabilities we can use the 



M{Cr) = -(0162 + &1&2 + aia2 - 6102) 



(39) 



«2 = j^i<^z 



oi = ax, 




(40) 
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formula [21] for the mean value of the index of criticality. However, now the mean value of 
this sum will be 2 a/2 (it is found straightforwardly) and for the mean value of the total 
number of non-critical impositions we obtain the value M{NonCr) = |(1 + ^ 0.85. 
The using of EPR pair of photons thus gives the substantial advantage for this problem. 

The effect we obtain could seem small, but if we repeat the assembly sequentially 
alternating it with some other processes, the probability can multiply and the effect may 
become very important. Moreover, this situation despite of its abstract character, in all 
likelihood is typical, and can appear in the different non-equilibrium processes in a living 
cell. 

For example, we look at the process of the penetration of a molecule through a cellular 
membrane. Two abstract factors a and b can determine the usefulness of it, which are the 
attributes of the cell itself and the penetrating molecule. For example, for the cell, it may 
be the different states of the molecule- helper, which pulls through the pore the considered 
molecule; for the penetrating molecule, it may be its chemical status. We associate with 
the cell the lower index 1, with the molecule the lower index 2. We assume that the sign 
determines is the pore open or not; for the molecule its state determining the readiness to 
the penetration through the pore. We then may build the construction analogous to the 
considered example, for which the bi-photonic control is more effective than the classical. 

We see that the application of the bi-photonic control in the different situations can 
raise the probability characteristic of efficiency comparatively with the classical control. 
Since the evolution of cell is the sequence of many steps, the probability characteristics 
of effectiveness will often multiply on all steps and the resulting effectiveness of the bi- 
photonic control can substantially exceed the classical analogue. 

The using of bi-photons for the establishing of the connections between the remote 
processors seems to be the most promising way of the development of the supercomputer 
technologies. This hybrid of classical silicon computer with the bi-photonic source is the 
most realistic embodiment of the idea of quantum computer nowadays. The development 
of this way can include also the more complex many photon states, for example, the 
entangled states of 4 and 6 photons that can help to represent more realistically the 
connection between processes simulated on the different computers. It is also possible 
to use the entangled states of particles with nonzero mass, for example, superconducting 
Cooper pairs, excited ion states in traps, or electrons in quantum dots. The software 
created on the basement of the collective behavior method makes possible to include 
easily entangled states of any type to our model. 

From the practical viewpoint, the using of entangled photon states is the best way 
to apply the miraculous quantum phenomenon of non-locality. I suppose that in the 
near perspective just this treatment of quantum computer is the most promising. The 
main area of applications for these computers will be the simulation of biochemistry. The 
further development of the idea of physical constructivism and its extension to the other 
types of complex systems is tightly connected with the development of the appropriate 
programming tools that requires the special efforts. 
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Figure 6: Constructive scheme of quantum computer 
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5 Conclusion 



As the conclusion of the course, I summarize the main features of physical constructivism. 

• Orientation to the dynamical scenarios instead of exact evaluations. 

• Pluralism: the admissibility of different scenarios, their comparison and selection. 
The direct including of the user into the model. 

• Possibility to consider mechanisms explicitly, on the level of samples. Scalability 
of the model. Preserving of the individuality of samples. Explicit form of the 
probability space. Explicit form of the entanglement as bonds between samples. 

The important component of the physical constructivism is the representation of dc- 
coherence as the fundamental factor, treated as the limitation of computational resources 
in the simulating computer. It brings the high price to the quantum entanglement effects 
obtained in experiments. Hence, the constructive approach to quantum computer differs 
from the Copenhagen (traditional) viewpoint. The quantum part of this computer is one 
of the controlling elements, e.g., is the kind of oracle. The constructive representation of 
quantum computer we schematically show at the picture. 

In the standard quantum theory, the heuristic is the intermediate step, which serves to 
the obtaining of equations, which solution is the main aim. In the constructive quantum 
theory, equations are the auxiliary tool for the obtaining of the heuristic, which then 
serves as the basis for the computer simulation. In the simulation of standard situations 
reducible to one particle, the constructive approach gives only the braking of modeling. 
It is effective only for the complex systems. What is the effect of constructivism for 
complex systems? It allows the direct manipulation with the entanglement as the physical 
phenomenon representing it as bonds. Standard Hilbert formalism has nothing similar. 
The presence of formal object bond opens new possibilities for us. With this formal 
apparatus, we hope to obtain the description of many atom chemistry, e.g., to advance 
on the way of investigation of complex systems. 

Somebody could object to me, that it would not be the physics at all. I would like to 
avoid the discussion about terms. The essence is only valuable: if we intend to investigate 
complex systems, we must pass in physics to the constructivism. Only in this case we could 
expect the real advance. The division of sciences follows from the presence of peculiar 
phenomena in each of them, and by advance, we mean the expansion of the exact methods 
to these phenomena, e.g., the expansion of physics. One should clearly understand that 
this expansion cannot be free of charge, and the cost of constructivism is the lowest one. 
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